Department of Optics 
Facultj? of Natural Sciences 

Palacky Universitj? 
Olomouc, Czech Republic 



Quantum properties of 
superposition states, 
squeezed states, and of some 
parametric processes 

Author: Faisal Aly Aly El-Orany 

Doctoral Thesis 
Branch: Optics and Optoelectronics 



Supervisor: Prof. RNDr. Jan Per ina, DrSc. 

Olomouc, 2001 



Acknowledgments 



First of all, I would like to start this thesis by passing my thanks to my supervisor. 
Professor Jan Pefina, under his supervision I worked more than three years; I thank 
him for collaboration in the most important and beautiful period of my life, I will never 
forget how much he was kind to me. I am very gratitude to him for long discussions 
and suggestions of interesting problems in nonlinear optics. This work could not be 
done in the present form without the support of all the members of the Department 
of Optics of Palacky University. 

My warm thanks come next to Professor Mohamed Sebawe Abdalla from the 
Mathematics Department, College of Science, King Saud University, Riyadh, Saudi 
Arabia, for his suggestions of several problems and helping in some cases in the 
calculations. I wish to thank him for everything he did for me from the earlier days 
when I started to be engaged with quantum optics. 

I would like to thank sincerely Professor Vlasta Perinova for many discussions, 
support, advise as well as supply with some articles. My acknowledge extends to Dr. 
Jaromir Kfepelka for many numerical advises. 

I am obliged much to my mother, brothers and sister, I missed them too much, 
although they are always with me, and to my wife and son who give beautiful meaning 
to my life, and support me and are patient with me. 



Doctoral Thesis 



2 



Faisal Aly Aly El-Orany 



C ontents 



1 Introduction 

2 Contemporary state of the problem 

2.1 Squeezed light 

2.1.1 Squeezed states and squeezed superposition states 

2.1.2 Entangled squeezed states 

2.1.3 Phase properties of squeezed states 

2.1.4 Squeezed states with thermal noise 

2.2 Light propagation 

2.2.1 Parametric processes 

2.2.2 Nonlinear directional coupler 

3 Goals of the thesis 

4 Methods and tools of quantum theory used 

4.1 Correlation functions 

4.2 Quadrature squeezing 

4.3 Quasiprobability functions 

4.4 Phase properties-Pegg-Barnett formalism 

5 Scientific results and their analysis: part I (static regime) 



1 



Contents 



5.1 Quantum statistical properties of superposition of squeezed and dis- 
placed number states 

5.1.1 State formalism and some of its properties 

5.1.2 Quasiprobability distribution functions 

5.1.3 Phase properties 

5.1.4 Scheme of production 

5.1.5 Conclusions 



31 



31 



37 



3S 



46 



45 



52 



5.2 Quantum statistical properties of superposition of squeezed and dis- 
placed states with thermal noise |51 

5.2.1 Density operator for the states |51 

5.2.2 Quasiprobability distribution functions 

5.2.3 Second-order correlation function |5E 

5.2.4 Photon-number distribution |56 

5.2.5 Phase property |5fi 

5.2.6 Conclusions 

5.3 Quantum statistics of a solvable three-boson squeeze operator model . 

5.3.1 Properties of the correlated quantum systems |6E 

5.3.2 Three-mode squeezed coherent and number states 

5.3.3 Quasiprobability distribution functions 

5.3.4 Conclusions |80 



6C 



63 



68 



74 



82 



6 Scientific results and their analjzsis: part II (dynamic regime) |81 

6.1 Quantum statistics and dynamics of nonlinear couplers with nonlinear 
exchange 

6.1.1 Squeezing phenomenon 

6.1.2 Quasiprobability functions |87 

6.1.3 Conclusions 

6.2 Generation of squeezed light in a nonlinear asymmetric directional coupler |93 



85 



91 



6.2.1 Equations of motion |93 



Doctoral Thesis 2 Faisal Aly Aly El-Orany 



Contents 



95 



6.2.2 Quasiprobability functions |96 

6.2.3 Photon-number distribution 

6.2.4 Phase distribution 

6.2.5 Conclusion 



100 



103 



6.3 Quantum statistical properties of nondegenerate optical parametric sym- 
metric coupler 

6.3.1 Equations of motion and their solutions 

6.3.2 Two-mode squeezing phenomenon 

6.3.3 Second-order correlation function 

6.3.4 Conclusions 



105 



106 



107 



109 



115 



7 Summary of main results 

8 Appendices 



117 



121 



Appendices Il20 
Bibliography 

Publications of the author of the thesis 1147 



131 



Strucne shrnuti vysledku disertacni prace (summary in Czech) Il49 



Doctoral Thesis 



3 



Faisal Aly Aly El-Orany 



Chapter 1 

Introduction 



The failure of classical mechanics to account for many experimental results such as the 
stability of atoms and matter, blackbody radiation, etc. led physicists to the realization of 
new tool to deal with these problems and thus quantum mechanics came to the vicinity 
of life at the beginning of the last century. Indeed, quantum mechanics is the one of the 
crowning achievements of modern physics where the quantized electromagnetic field 
is supported by the experimental observations of nonclassical states of the radiation 
field, e.g. squeezed states, sub-Poissonian photon statistics and photon antibunching. 
Nowadays, quantum optics, the union of quantum field theory and physical optics, is 
undergoing a time of revolutionary change. The subject has evolved from early studies 
on the coherence properties of radiation to the laser in the modern areas of study 
involving, e.g., the role of squeezed states of the radiation field and atomic coherence 
in quenching quantum noise in interferometry and optical amplifiers. Furthermore, 
quantum optics provides a powerful new probe for addressing fundamental issues 
of quantum mechanics such as complementarity, hidden variables, and other aspects 
central to the foundations of quantum physics and philosophy. 

For the description of a quantum system, the concept of state is used, which is 
the same, a wave function, a state vector or a density matrix containing the inform- 
ation about the possible results of measurement on the system. Quantum optics has 
statistical origin and therefore the state of a quantum system contains all information 
necessary to completely determine its statistics (the probabilistic nature of a quantum 
system). Among the wide variety of possible radiation-field states, there are some 
fundamental types that play a special role in quantum optics. Strictly speaking, there 
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are three types of these states which are widely used, namely, coherent, number and 
squeezed states. Following the development of the quantum theory of radiation and 
with the advent of the laser, the coherent states of the field, that mostly describe a 
classical electromagnetic field, were widely studied. Indeed, these states are more 
appropriate basis for many optical fields. However, number states are purely nonclas- 
sical states and they are a useful representation of high-energy photons. However, 
there are experimental difficulties which have prevented the generation of photon 
number states with more than a small number of photons. Despite this the number 
states of the electromagnetic field have been used as a basis for several problems in 
quantum optics including some laser theories. The squeezed states (i.e. the states of 
light with reduced fluctuations in one quadrature below the level associated with the 
vacuum state) are very important owing to their potential applications, e.g. in the op- 
tical communication systems, interferometric techniques, and in an optical waveguide 
tap. Furthermore, squeezed states have been seen in several laboratories. In addition 
to these states, several quantum states in the literature have appeared, in particular 
those which are connected with the superposition principle. Using such principle, in 
the first part of this thesis, we will develop new class of states which are superposition 
of squeezed and displaced number states. In fact, this class of states generalizes some 
considerable results given in the literature. We close this part by giving an example 
of multidimensional squeeze operator which is more general than usually used. 

On the other hand, nonlinear optics (NLO) has become a very important subfield of 
optics since its inception over 40 years ago. The origin of this branch is the study of the 
phenomena that occur as a consequence of the modification of the optical properties 
of a material system by the presence of light. The impact of NLO on science and 
technology has been twofold. First, it has enhanced our understanding of fundamental 
light-matter interactions. Second, it has been a deriving force in the rejuvenation of 
optical technology for several areas of science and engineering. For example, NLO 
provides the key to many developments, e.g. the advent of telecommunications using 
optical fibers, carrying information on a laser beam in the process of communication, 
storage, retrieval, printing or sensing, and there are increasing efforts to achieve ever 
greater data-processing capabilities. One of the most promising devices in NLO is the 
nonlinear directional coupler that consists of two or more parallel optical waveguides 
fabricated from some nonlinear material. Both waveguides are placed close enough 
to permit flux-dependent transfer of energy between them. This flux transfer can be 
controlled by the device design and the input flux as well. This device is experimentally 
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implemented. Several models of this device will be investigated in the second part of 
this thesis. 

As is well-known the quantum statistical properties of radiation present an import- 
ant branch of modern physics with rapidly increasing applications in spectroscopy, 
quantum generators of radiation, optical communication, etc. Using these methods, 
quantum optics was able to predict many nonclassical phenomena, such as squeez- 
ing of vacuum fluctuations in field quadratures, antibunching of photons, sub-Poisson 
photon statistics exhibiting reduced photon-number fluctuations below the Poisson 
level of ideal laser light, quantum oscillations in photocount distributions, violation of 
various classical statistical inequalities, collapse and revival of atomic inversion, etc. 

This doctoral thesis deals with investigating various quantum statistical aspects for 
new quantum states and for nonlinear directional couplers. The first part will be 
devoted to the static regime where we introduce the superposition of squeezed and 
displaced number states and discuss their quantum properties, phase properties and 
the influence of thermal noise on their behaviour. We conclude this part by giving 
an example of new type of multidimensional squeeze operator model which is more 
general than usually used and which includes two different squeezing mechanisms. 
The dynamical regime will be considered in the second part where we study the evol- 
ution of quantum states like coherent, number and thermal states in some nonlinear 
medium, such as nonlinear couplers. To achieve this goal in a systematic way we give 
first a survey of previous works in the literature related to our topic in chapter 2 and in 
chapter 3 we write down all mathematical relations controlling the quantum statistical 
properties of the models under discussion. 
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Chapter 2 

Contemporarj; state of the problem 



In this chapter, we display the most important relevant results in quantum optics which 
are used in our work. This will be done by throwing the light on the results of both 
squeezed light and the propagation of light in the nonlinear media. 

2.1 Squeezed light 

2.1.1 Squeezed states and squeezed superposition states 

Squeezed states of the electromagnetic field are purely quantum states. They are 
defined through the relation 

|ei,a) = D(a)S(ei)|0), (2.1.1) 

where |0) is the vacuum state, D(a) [l] and S(^i) |^ are displacement and squeeze 
operators which are given, respectively, by 

D(a) = exp(d''a - da*), (2.1.2) 

S(ei) = exp(ei*d2 - ^id^), (2.1.3) 

where d and d ^ are annihilation and creation operators, respectively, and and a are 
complex parameters. 



9 



Chapter 2. Contemporarv state of the problem 



The significant feature for these states is that they can have less uncertainty in 
one quadrature than a coherent state. These states can exhibit a number of distinctly 
quantum features, such as sub-Poissonian statistics [3]-l5], as well as they have no 
nonsingular representation in terms of the Glauber-Sudarshan P distribution [6]. It is 
important to refer that the concept of squeezed coherent states (SCS) have been applied 
to other quantum mechanical systems. For instance, they may play a role in increasing 
the sensitivity of a gravitational wave detector SCS have appeared at the first time 
[2] as a simple generalization of the well known minimum-uncertainty wavepackets. 
Authors of |[6][7] have demonstrated further details concerning the properties of such 
states, or the so-called two-photon coherent states [7]. 

A lot of efforts have been directed towards the methods of generating SCS. For 
example, we can mention, authors of P-[TOl have shown that squeezing may be gen- 
erated in an optical four-wave mixing. Further progress has been done in papers 
llT][l2l, where the first experimental observations of squeezing were given in four- 
wave mixing in an atomic beam of Sodium vapor (11] as well as in optical fibers 
(12) . Particular attention was given to generate squeezed light in parametric amplifiers 
P ITStiTB] . where the treatment of parametric oscillation and intracavity second har- 
monic generation provides the basis for subsequent calculations of squeezing in these 
systems. Furthermore, more activities have been focused on studying generation of 
SCS in: nondegenerate parametric oscillator [17], optical bistability Il8l[l9], and res- 
onance fluorescence j20l - [22] . For more complete information about SCS the reader 
can consult papers p3l - [25| . 

The concept of superposition of quantum states has been extended by several au- 
thors to include SCS [261127] . Authors of [26] have introduced even and odd displaced 
squeezed states as a superposition of two vacuum squeezed displaced states. For such 
superposition the authors studied the higher-order squeezing with respect to the defin- 
ition given by Hong and Mandel [28]. Also it is interesting to refer to ^3 where 
extensive efforts have been done to introduce a superposition of set of SCS and the 
authors managed to calculate and discuss their quasiprobabality distribution functions 
as well as the generation scheme of such superposition states. 

Squeezed and displaced number state (SDNS) [29] is an energy eigenstate of a 
quantum harmonic oscillator which is displaced and then squeezed (it may also be 
called the generalized squeezed coherent state (GSCS)). For SDNS bunching and anti- 
bunching properties have been investigated ^Ol- For more details about its properties, 
discussion as well as new methods for analytical investigation one can see paper [51] . 



Doctoral Thesis 



10 



Faisal Aly Aly El-Orany 



2.1. Squeezed light 



Further, for such a state, the most general form as well as the time-dependent ex- 
pectation values, uncertainties of wave-functions and probability densities have been 
given in ||32] using the functional form for the squeeze and the time-displacement 
operators j33]. Also we can mention that a special case of SDNS has been discussed 
in |34] where these states have been treated from the point of view of non-diagonal 
P-representation. The authors have discussed the properties of the squeezed and dis- 
placed Fock states as generalized states and their discussion has been extended to the 
Glauber R-representation for the density operator as well as to the phase distribution. 
The physical interpretation of SDNS has been considered similarly to SCS |29], in 
other words, SDNS is the coherent state formed due to two excitations on a particular 
number state. Also we can mention paper |35] where squeezed (but not displaced) 
number state was produced. 

We can also refer to several applications of SCS. For example, SCS have several 
potential applications, one, for instance, is in optical communication systems 136]. Also 
the interferometric techniques may provide ways to detect very weak forces such as 
caused by gravitational radiation and may experience limitations on sensitivity due to 
quantum noise arising from photon counting and radiation pressure fluctuations |37J. 
Another application is in an optical waveguide tap where it has been shown that 
a high signal-to-noise ratio may be obtained using SCS in an optical waveguide to tap 
a signal carrying waveguide. 



2.1.2 Entangled squeezed states 

Entangled states gain their feature from the quantum correlation between different 
quantum mechanical systems where the individual operator of the single system does 
not exhibit nonclassical effects, however, the compound system can exhibit such effects. 
To be more specific, if the measurement of an observable of the first system (say), for 
correlated system, is performed, this projects the other system into new states; other- 
wise the systems are uncorrected. Squeezing property is the important phenomenon 
distinguishing well mechanism of correlation of systems, where squeezing can occur 
in combination of the quantum mechanical systems even if single systems are not 
themselves squeezed. 

The most significant example for compound squeezing is related to the two-mode 
squeeze operator defined by an effective unitary operator as 
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5(^2) = exp(e2*d^b^ - ^206), (2.1.4) 

where d and b are the annihilation operators of first and second modes respectively, 
and ^2 is a complex parameters. Expression ( |2.1.4| represent the evolution operator 
of the nondegenerate parametric amplifier with classical pumping and this operator 
can produce perfect squeezing only in the correlated states of two field modes. There 
are several models which have been given in the literature, e.g. I39]- |43]. The aim 
of these articles is to examine the statistical properties for single- and compound- 
modes. Moreover, the ideas that quantum correlation can give rise to squeezing in the 
combination of system operators has been shown true for multimode squeezed states 
of light |4I5II11I39II44II45] and for dipole fluctuations in multimode squeezed states f40] . 

2.1.3 Phase properties of squeezed states 

In classical optics, the concepts of the intensity and phase of optical fields have a 
well-defined meaning. That is the electromagnetic field (E) associated with one mode, 
E = Aexp(i9), has a well defined amplitude (A) and phase (9). This is not so simple 
in quantum optics where the mean photon number and the phase are represented 
by noncommuting operators and consequently they cannot be defined well simultan- 
eously. In fact, the concept of phase is a controversional problem from the earlier days 
of quantum optics (461 - 148] . In general there are three methods of treating this issue 
149]. The first one considers the phase as a physical quantity in analogy to position 
or momentum by representing it with a linear Hermitian phase operator. The second 
one involves c-number variables (real or complex) in phase spaces or their associated 
distribution functions, or ensembles of trajectories. The third one is the operational 
phase approach in which the phase information is inferred from the experimental data 
by analogy with the classical analysis of the experiment. Each approach has advantage 
and disadvantage points. 

Squeezed states have phase sensitive noise properties and therefore several works 
have been devoted to follow such properties. We can mention that the authors of |j50] 
have investigated the fluctuation properties of squeezed states using a phase-operator 
formalism defined by Susskind and Glogower [51] • They have shown that similarly as 
coherent states of high intensity approach a semiclassical number-phase uncertainty 
product, the squeezed states retain their quantum properties and their number-phase 
uncertainty relations are not minimized. However, the exact general results of the 
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same technique have been obtained in f52][53] for single- and two-mode squeezed 
states. Indeed, these are mathematical treatment for the problem. Exact phase calcu- 
lations for different definitions have been presented in ^4] showing that the measured- 
phase-operator formalism leads to contrasting behaviour compared with that based on 
the Susskind-Glogower or Hermitian-phase-operator formalisms. In the framework 
of Pegg-Barnett formalism f55] - l58] several works have been done treating not only 
single-mode squeezed states f59ll60] but also two-mode squeezed states |61J- The main 
results of these articles are: for single-mode squeezed states with non-zero displace- 
ment coherent amplitude, the phase distribution exhibits the bifurcation phenomenon 
or single peak-structure under certain conditions. However, it has been shown that 
the joint phase distribution for the two-mode squeezed vacuum depends only on the 
sum of the phases of the two modes, and that the sum of the two phases is locked to 
a certain value as the squeeze parameter increases f62] . 

We refer to j47l48] where the investigation of phase properties of various quantum 
mechanical models in greater details have been given. 

2.1.4 Squeezed states with thermal noise 

It is worthwhile mentioning that squeezing of thermal radiation field has been already 
produced in a microwave Josephson-junction parametric amplifier |j63] where a thermal 
input field has been considered in a squeezing device and the generated field exhibits 
substantial noise reduction. The aim of such work is to generate nonclassical fields for 
interaction studies with Rydberg atoms in high-Q microwave cavities, where thermal 
noise in input fields is always large. Further, it has been suggested [37] that inter- 
ferometers for detection of gravitational waves could employ squeezing techniques in 
order to improve their resolution. Practically, these systems will inevitably experience 
thermal noise, so it is essential to be aware of the various representations of squeezed 
states and their physical interpretation j64]. 

As a result of the fact that signal beams are usually accompanied by thermal noise, 
many authors concentrated on the studying of the influence of thermal noise on the 
behaviour of quantum states |65]- |[75]. Some of these studies give particular attention 
to the calculation of photon-counting distribution [65]- |[69]. This calculation basically 
depends on the Gaussian form of Wigner function. It has been shown that for large 
values of the squeeze parameter the photon-counting distribution is oscillating even 
for strong thermal noise. The quantum statistical properties of squeezed thermal light 
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f70)- | [74] , e.g. quadrature squeezing, second- and higher-order correlation functions, 
different representation for the density operator, and quasiprobability distribution func- 
tions, reveal that the degree of purity of the input thermal light is left unchanged by the 
subsequent squeezing and displacement processes. Moreover, an effective squeezing 
could be recognized under particular choice of the parameters and this fact is ensured 
by the behaviour of Glauber P-function which it no longer exists. For strong squeez- 
ing both squeezing properties and normalized correlation functions to all orders do 
not depend on the initial intensity of thermal field. The quantum phase distributions 
and variances of strong coherent and phase-squeezed states mixed with thermal light 
have been considered in ||75] showing that the effect of thermal noise on the coherent 
phase distribution becomes important only when the number of thermal photons is of 
the order of one-half ((l/2)e"'', r is the real squeeze parameter). 

Here we may refer also to the squeezed thermal reservoir which has been studied 
in detail 126]- IZS]- One can see that the authors of ref. [78] have shown for atom 
radiation in a squeezed thermal reservoir that the two quadratures of the atomic po- 
larization are damped at different rates, which is consistent with the case of a squeezed 
vacuum reservoir. Furthermore, when the atom is driven by a coherent field, it has 
been found that the steady-state polarization depends on the relative phase of the 
squeezing and driving field. This phase dependence becomes less pronounced when 
the number of thermal photons increases. This behaviour suggests a number of novel 
applications, such as new schemes for optical bistability |80]. Nevertheless, an objec- 
tion against the reservoir technique can be given when an exact solution of the density 
matrix equation becomes unavailable even in the steady-state regime. 



2.2 Light propagation 

Classical optics describes quite successfully the propagation of laser light, both in 
free-space and inside a transparent medium. The reason is that in a coherent state of 
radiation, the electric and magnetic fields may be written in terms of their expectation 
values, and thus their propagation may be treated classically through the macroscopic 
Maxwell equations. The classical Maxwell equations permit the calculation of both 
the spatial progression and the temporal evolution of a propagating electromagnetic 
field, and treat the interaction of the field in a medium phenomenologically through 
the induced polarization. However, in quantum optics the simultaneous use of the 
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Hamiltonian and the momentum operators yields operatorial spatial-temporal equa- 
tions of motion for the electric and magnetic fields, having a form equivalent to that 
of the classical Maxwell equations and thus the quantum propagative phenomena can 
be rigorously described [811182] , 

In this section we throw the light on two examples representing the propagation 
of light in media and which are frequently used in this thesis; namely, parametric 
processes and nonlinear directional coupler. 

2.2.1 Parametric processes 

There has been a great interest in the field of nonlinear optics for both the practical 
applications and the theoretical aspects of the nonlinear effects. Experiments in this 
field were made possible by the fact that the lasers with a sufficiently high output (10^ 
to 10^ Watt/cm^) had become easily available [83]. At this power level the nonlinear 
susceptibilities of certain media were producing observable effects |j84], e.g. such as 
the phenomena of parametric fluorescence and parametric oscillation ||85]. Indeed, 
linearity or nonlinearity is a property of the medium through which the light travels, 
rather than the property of the light itself. Light interacts with light via the medium. 
More precisely, the presence of an optical field modifies the properties of the medium 
which, in turn, modify another optical field or even the original field itself pS]. For 
example, the study of the structure of crystals and molecules has long utilized the phe- 
nomena of light scattering from atoms or molecules having two energy levels. The 
frequency of the incident beam may then be shifted up or down, by an amount equal 
to the difference in the two energy levels of the scatterer. The resulting lower- and 
higher- frequency scattered waves are the Stokes and anti-Stokes components, respect- 
ively |87j. However, in the coherent Raman effect the presence of a monochromatic 
light wave in a Raman active medium gives rise to parametric coupling between an op- 
tical vibrational mode and a mode of the radiation field which represents the scattered 
(Stokes) wave. In the case of Brillouin scattering a similar form of coupling holds, with 
the vibrational mode oscillating at an acoustic rather than an optical frequency |5]. 

As was known the nonlinear processes in the quantum mechanical domain has led 
to the prediction and the observation of many quantum phenomena, e.g. squeezing of 
vacuum fluctuations and photon antibunching. In the heart of nonlinear optics there 
are two significant processes which have been attracted amount of study; namely, 
parametric frequency converter (PFC) and parametric amplifier (PA). PEC can be de- 
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scribed by a process of exchanging photons between two optical fields of different 
frequencies: signal mode at frequency and idler mode at frequency (i»2- This model 
can be applied to describe various optical phenomena, e.g. to find analogies between 
PFC and beam splitter |^8], two-level atom driven by a single mode of electromagnetic 
field |j89], and Raman scattering f88l90] . The quantum properties of PFC are discussed 
in j9l]. Further, some authors studied this model as the lossless linear coupler, e.g. 
(92] - (95] . In this situation the model is considered to be represented by two electro- 
magnetic waves which are guided inside the structure consisting of two adjacent and 
parallel waveguides; the linear exchange of energy between these two waveguides is 
established via the evanescent field | j96] . 

On the other hand, PA is designed in the most familiar form to amplify an oscillating 
signal by means of a particular coupling of the mode to a second mode of oscillation, 
the idler mode. The coupling parameter is made to oscillate with time in a way which 
gives rise to a steady increase of the energy in both the signal and idler modes j5]. 
The importance of PA is related to the fact that it is the source for squeezed light (7]. 
For example, degenerate and non-degenerate PA are sources for single-mode fT] and 
two-mode f39] - l4T] squeezing of vacuum fluctuations, respectively. 

The parametric process has been employed in experiments. For example, the 
fourth-order interference effects arise when pairs of photons produced in paramet- 
ric down-conversion are injected into Michelson interferometers |[97]. However, the 
second-order interference is observed in the superposition of signal photons from two 
coherently pumped parametric down-conversions when the paths of the idler photons 
are aligned pSj. Further, squeezed states of the electromagnetic field are generated 
by degenerate parametric down-conversion in optical cavity [16] where noise reduc- 
tions greater than 50% relative to the vacuum noise level are observed in a balanced 
homodyne detector. Also, the observation of high-intensity sub-Poissonian light using 
the correlated "twin" beams generated by an optical parametric oscillator has been 
demonstrated (99] . 

The optical processes involving the competition between PFC and PA are of interest 
from theoretical and experimental points of view, e.g. in three-mode interaction (87) 
1100] . For these processes the quantum theory can be briefly reported as follows (101] : 
The nonlinear interaction couples different photon modes and leads to energy transfer 
between modes. Photons in some modes may be annihilated, while those in the other 
modes created, and hence the photon distribution is disturbed. In every time the rate 
of energy transfer between the modes depends on the statistical properties of the 
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light fields. Statistics is particularly important in this case for analyzing the results of 
experiments. As is expected the output field depends on the statistical nature of both 
the incident beams and the fluctuations in the medium. Moreover, the measurements 
of the statistical properties of the output field could yield information on the properties 
of the medium if those of the incident radiation are known. 

2.2.2 Nonlinear directional coupler 

In quantum optics many simple quantum systems have been examined from the point 
of view of completely quantum statistical description including not only amplitude and 
intensity (energy) development of such systems, but also higher-order moments and 
complete statistical behaviour. Such results have fundamental physical meaning for 
interpretation of quantum theory jl02] and they are useful for applications in opto- 
electronics and photonics as well. These results can be successfully transferred to 
more complicated and more practical systems, such as optical couplers composed of 
two or more waveguides connected linearly by means of evanescent waves. The wave- 
guides can be linear or nonlinear employing various nonlinear optical processes, such 
as optical parametric processes, Kerr effect, Raman or Brillouin scattering, etc. Such 
devices play important role in optics, optoelectronics and photonics as switching and 
memory elements for all-optical devices (optical processors and computers). When 
one linear and the other nonlinear waveguides are employed, we have a nonlinear 
optical coupler producing nonclassical light in the nonlinear waveguide which can be 
controlled from the linear waveguide, i.e. one can control light by light. The gener- 
ation and transmission of nonclassical light exhibiting squeezed vacuum fluctuations 
and/or sub-Poissonian photon statistics in nonlinear optical couplers can further be 
supported when all the waveguides are nonlinear. The possibility to generate and 
to transmit effectively nonclassical light in this way is interesting especially in optical 
communication and high-precision measurements where the reduction of quantum 
noise increases the precision. 

Since the pioneering work on nonlinear couplers which has been done by Jensen 
[105] , a series of articles have been devoted to the study of this important optical device 
from both classical fl04j - [114] (using the coupled-mode theory) and quantal viewpoints 
[115] - [155] . In the framework of quantum mechanics particular attention has been 
paid to quantum statistical properties in relation to quantum noise and generation and 
transmission of nonclassical light. Such generation and transmission of nonclassical 
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light can be very effective as a consequence of using evanescent waves involved in 
the interaction. For a review of role of quantum statistical properties in nonlinear 
couplers, see [136] . 
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Goals of the thesis 



The main goal of the doctoral thesis is a study of quantum statistics for some static 
and dynamic regimes of nonlinear processes in quantum optics. 

In the first part (static regime): 1- we want to develop a general class of quantum 
states as a superposition of displaced and squeezed number states. We study the quan- 
tum statistics for this class rigorously. We also investigate the effect of thermal noise 
on the properties of such class of states. Also we report the methods of generation 
for such superposition. 

2- we want to give an example of new type of multidimensional squeeze operator model 
which is more general than usually used and which includes two different squeezing 
mechanisms. All basic properties related to this operator are discussed in greater 
details. 

In the second part of this thesis we concentrate on studying the statistical properties 
of an optical field propagating within a nonlinear directional coupler. Our starting 
point is the Hamiltonian, which represents a nonlinear directional coupler. We assume 
one-passage propagation, so that losses in the beams have been neglected. In other 
cases they can be described in the standard quantum way in the form of interaction 
of light beams with reservoirs, as for instance described in |157] . Moreover, we 
treat the problems of propagation in the Hamiltonian formalism assuming the energy 
of the system does not have directionality. However, in the case that all waves are 
propagating with the same velocity, time and space relate by the velocity of propagation 
V, z = vf. For all these models we investigate the effect of switching between the input 
modes and the outgoing fields from the coupler. We can consider three problems 
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of light propagation in this device: 1- a symmetric directional coupler operating by 
nondegenerate parametric amplification. 

2- nonlinear optical couplers composed of two nonlinear waveguides operating by the 
second subharmonic generation, which are coupled linearly through evanescent waves 
and nonlinearly through nondegenerate optical parametric interaction. 

3- a nonlinear asymmetric directional coupler composed of a linear waveguide and a 
nonlinear waveguide operating by nondegenerate parametric amplification. 
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Chapter 4 

Methods and tools of quantum 
theorj; used 

In this chapter we review the quantum methods and tools for controlling the nonclas- 
sical phenomena, such as correlation function, quadratures squeezing, quasiprobability 
distribution function, photon-number distribution and phase distribution. Further, we 
do not describe a complete details for these methods and tools when they can be 
found in the standard text book. Moreover, we write down the formulae related to the 
dynamical regime where those of the static regime can be obtained by simply setting 
f = as a main tool the coherent-state technique is used. 

4.1 Correlation functions 

Antibunched and/or sub-Poissonian light is an example of nonclassical light and can 
be determined from a photocounting-correlation measurement. Starting with the ex- 
periment of Hanbury Brown and Twiss |158] , strong interest in the photon-counting 
statistics of optical fields began. Traditional diffraction and interference experiments 
and spectral measurements may be considered as being performed in the domain 
of one photon or linear optics. The theory of higher-order optical phenomena, de- 
scribed by higher-order correlation functions of the electromagnetic field, was founded 
by Glauber |139| , who introduced the measure of super-Poissonian statistics (classical 
phenomenon) and sub-Poissonian statistics (nonclassical phenomenon) of photons in 
any state. A state (of a single mode for convenience) which displays sub-Poisson stat- 
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istics is characterized by the fact that the variance of the photon number {{Ahj{t))^) 
is less than the average photon number (n;(f)) = {Aj{'t)Aj{t)). This can be written by 
means of the normalized normal second-order correlation function as 



(2j (af(f)d2(f)) ((An,-(f))2) - (df(f)d,-(f)) 

g. '(t) = 'rr = 1 + X , ('i-l-l 

' (d/(f)d,-(f))2 (d/(f)d,-(f))2 

where the subscript ; relates to the ;th mode, and the photon number variances have 
the form 



((An,-(f))2) = ((d/(f)d,-(f))2) - (d/(f)d,-(f))2. (4.1.2) 

(9) (9) 

Then it holds that g- (f) < 1 for sub-Poissonian distribution of photons, gy (f) > 1 for 
super-Poissonian distribution of photons and when (f) = 1 Poissonian distribution 
occurs corresponding to a coherent state. Furthermore, for instance, the generation 
of sub-Poissonian light has been established in a semiconductor laser |140j and in the 
microwave region using masers operating in the microscopic regime jl41] . An ap- 
plication of radiation exhibiting the sub-Poissonian statistics to optical communications 
has been considered in |142] . 

On the other hand, it has been shown explicitly in fl43lll44] that sub-Poissonian 
photon statistics need not be associated with antibunching, but can be accompanied 
by bunching. However, within the framework of the classical theory, light cannot be 
antibunched, i.e. antibunched light is a manifestation of a quantum effect. The basic 
formula to study this phenomenon is the two-time normalized intensity correlation 
function |144II145] . For the 7th mode, this function is defined by 

, ^ ,) . <f f ^ ^"^'"» . ,4.1.3) 

' (d|(f)d,-(f))(d|(f -h T)d,-(f -h t)) 

The importance of this function in the analysis of photon antibunching comes from the 
direct relation between this function and the joint detection probability of two photons, 
one at time t and another at time f + t. It is clear that using l |4.1.3| for t ^ as 
a definition of bunching properties, then the bunching/antibunching and super-/sub- 
Poissonian statistics are in one-to-one correspondence. More general definition of 
photon antibunching can be adopted |i 44111 45| if gf\t, i + t) increases from its initial 
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value at T = 0. This can be represented in equivalent differential form, assuming that 

(9) 

Qj [t.t + t) is a well behaved function in t, as 

dgf\t,t + t) 

Ht) = -^^-^ ^|r=o > 0; (4.1.4) 

photon bunching is given by the opposite condition (K/(f) < 0), otherwise the photons 
are unbunched. It is reasonable pointing out that antibunching and sub-Poissonian 
behaviour always accompany each other for the single-mode time-independent fields 

mi- 

Finally, we turn our attention to discuss the effect of intermodal correlation in 
terms of anticorrelations between different modes in the model. This can be done 
by two means. The first mean is given by introducing the photon-number operator 
nyjj = A|Ay + A^Ajj and calculating the quantity ((AW,jj)^) = 

(: [nj.kf '■) - {fi-j,k)^> where : : denotes the normally ordered operator, i.e. creation 
operators A| are to the left of annihilation operators Ay. The quantum anticorrelation 
effect is then characterized in terms of the variance of the photon number, which is 
less than the average of the photon number for nonclassical light, by negative values of 
((A W,jj)^), i.e. negative cross-correlation taken two times is stronger than the sum of 
quantum noise in single modes jl47j . The second way is based on violation of Cauchy- 
Schwarz inequality. The violation of Cauchy-Schwarz inequality can be observed in a 
two-photon interference experiment |148] . Classically, Cauchy-Schwarz inequality has 
the form (1491 



{hhf < (7f)(7|), (4.1.5) 

where Jy,; = 1,2 are classical intensities of light measured by different detectors in a 
double-beam experiment. In quantum theory, the deviation from this classical inequal- 
ity can be represented by the factor |150| 



[(A["A2)(AfA|)]i 
(A|AyA[A,) 

The negative values for the quantity Ij ^ mean that the intermodal correlation is larger 
than the correlation between the photons in the same mode f43] and this indicates 
strong violation of the Cauchy-Schwarz inequality. Finally, anticorrelation between 
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modes can be measured by detecting single modes separately by two photodetectors 
and correlating their outputs. 

4.2 Quadrature squeezing 

As we mentioned earlier, squeezed light possesses less noise than a coherent light in 
one of the field quadratures and can exhibit a number of features having no classical 
analogue. This light can be measured by homodyne detection where the signal is 
superimposed on a strong coherent beam of the local oscillator. 

There are several definitions for squeezing, e.g. standard squeezing jl51] , amplitude- 
squared squeezing |152] , two-mode squeezing |155] , higher-order squeezing |154] , 
principal squeezing |155] , etc. Of course, the quantum mechanical systems can ex- 
hibit different types of squeezing at the same time. It has been shown that when a 
beam of light propagates through a nonlinear crystal, in the process of generation 
of the second harmonic, the fundamental mode becomes squeezed in the sense of 
standard [151j as well as amplitude-squared squeezing fl52] . Further, the parametric 
amplifier in a cavity can ideally produce squeezed light with characteristics akin to 
the single and two modes when operating in degenerate and non-degenerate regimes, 
respectively [155] . 

In this thesis we investigate single-mode, two-mode and three-mode squeezing on 
the basis of the two quadratures and (where the subscript n takes on 1,2,3 
associated with the single-, two- and three-mode squeezing) which are related to the 
conjugate electric and magnetic field operators E and H. They are defined in the 
standard way. Assuming that these two quadrature operators satisfy the following 
commutation relation 



where C is a c-number specified later, the following uncertainty relation holds 



where {{AXnf) = (X^) _ (x^)2 the variance. Therefore, we can say that the model 
possesses Xn-quadrature squeezing if the S^-factor [156] , 



(4.2.1) 



((AX„)2)((AV„)2) > 




(4.2.2) 
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Sn = ((AX„(f))2) -0.5|C|, 
((AX„(f))2)-0.5|C| 



(4.2.3) 



0.5|C| 

satisfies the inequality -1 < < 0. Similar expression for the -quadrature (Q^- 
parameter) can be obtained. 

For example, the two quadratures of three-mode squeezing are defined as 

X^{t) = ^[Ai(f) + A2(f) + A3(f) + A[(f) + Alit) + Alit)], (4.2.4) 

= ^[Ai(f) + A2(f) + A3(f) - A[(f) - Alit) - Alit)]. (4.2.5) 

The expressions for the single-mode and two-mode squeezing can be obtained easily 
from j4.2.4| and l |4.2.5| by dropping the operators of absent mode, e.g. for the 1st 
mode, single-mode squeezing can be obtained by setting the operators of 2nd and 3rd 
modes equal zero. It should be taken into account that C = 5/ 1/ 1 corresponding to 
the single-mode, two-mode and three-mode squeezing, respectively. 



4.3 Quasiprobabilifj; functions 



Evaluation of various time-dependent mode observables is most conveniently achieved 
with the aid of corresponding time-dependent characteristic functions, their normal, 
antinormal and symmetric forms, and the Fourier transforms of these characteristic 
functions (quasiprobability functions). All of these are related to the density matrix 
which provides a complete statistical description of the system. There are three types 
of quasiprobability functions: Wigner W-, Glauber P-, and Husimi Q-functions. These 
functions could be used also as crucial to describe the nonclassical effects of the system, 
e.g. one can employ the negative values of W-function, stretching of Q-function and 
high singularities in P-function. Furthermore, these functions are now accessible from 
measurements |157] . 

Indeed, the detailed statistics of the coupled field modes can be obtained from 
several photon-counting experiments. Most often we are interested in the quantum 
statistics of either one mode which determine the ensemble averages of the observable 
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of this mode, or the composite statistics of the compound modes which reflect their 
mutually correlated properties. That is why we consider in this thesis phase space 
distributions for the single- and compound-modes for different types of initial states. 

The s-parametrized joint characteristic function for the compound system (n modes) 
defined by 



d")(C^s) = Tr 



p(0) exp 



^(oAf(f)-c;A;(f) + -ioi 

7 = 1 



(4.3.1) 



where p(0) is the initial density matrix for the system under consideration, C = (Ci- C2' Cn) 
is a set of parameters, and s takes on values 1,0 and -1 corresponding to normally, 
symmetrically and antinormally ordered characteristic functions, respectively. 

Thus the s-parametrized joint quasiprobability functions are given by 

W^^^ifi'ts) = -L [ C(")(CJ,s)f]exp(j8,C; -J8;0)d'0- (4.3.2) 

— 7 = 1 

where j8 = {Pi,P2> ■••>Pn) are complex amplitudes and / = / .../ n-fold integral. When 
s = 1,0,-1, equation (|4.3.2t gives formally Glauber P-function, Wigner W-function 
and Husimi Q-function, respectively. 

On the other hand, the s-parametrized quasiprobability functions for the single- 
mode can be attained by means of integrating (n - 1) times in the corresponding joint 
quasiprobability functions or by using the single-mode s-parametrized characteristic 
function. For instance, the characteristic function C'^)(Ci,f,s) for mode A\ (say) can 
be obtained from C'")(C, ^,s) by simply setting C,k = 0, k = 2, 3, ...,n. Hence the 
s-parametrized single-mode quasiprobability can be derived by 



r " 
wW(j8i,f,s)= W^''^{l,t,s)Y]dpj, 

— 7 = 2 



(4.3.3) 



etc., or 



wW(i8;,f,s) = I C^'HCi,t,s)exp{PjCj -PjCi)d% 



(4.3.4) 



C^^\C,j, t, s) in ( |4.3.4| is the s-parametrized characteristic function for single-mode. The 



Doctoral Thesis 



26 



Faisal Aly Aly El-Orany 



4.3. Quasiprohahility functions 



superscripts (1) and (n) in the above equations stand for single-mode case and n-mode 
case, respectively. 

The various moments of the bosonic operators for the system, using the charac- 
teristic functions and quasiprobability functions, in the normal form (N), antinormal 
form (A) and symmetrical form (S), corresponding to s = 1,-1,0, respectively, can be 
obtained by 

fi 

where n,-, my are positive integers. The formulae l |4.5.5t are valid for the single mode 
and compound modes as well. 

There are several applications of the quasiprobability function. We restrict ourselves 
to those of the single-mode case because the generalization to the multi-mode case is 
straightforward. 

First, the degree of purity Trp^ for any state can be evaluated with aid of the 
symmetrical characteristic function given by the relation 

Tr(p2) = l I |C(C,s = 0)|2d2c. (4.3.6) 

One can find for a pure state Tr(p^) = 1 and for a mixed state Tr(p^) < 1. Also, the 
quasiprobability function can be used to calculate the matrix elements of pi(f) in the 
n-quantum representation at time t by means of the relation |5] 

/nmn*ni 
Pl^w'^'^ifi^.s = l,f)exp(-|i8i|Vi8i. (4.3.7) 

It is clear that this formula reduces to the photon-number distribution when m = n\. 
Further, the photon-number distribution can be obtained also from the integral relation 
in terms of Wigner function of the specified mode and Laguerre polynomials as 

P^'\n^.t) = I W(i)(i8,-,f)exp(-2|i8,f )L„.(4|i8,f )d2j8,-, (4.3.8) 
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where ; denotes the mode under consideration and W'^^(j8y,f) is the single mode 
Wigner function for the specified mode. Also quasiprobability functions can be used 
to investigate the phase distribution for the structure by integrating these functions 
over the radial variable 



roo 

P{e,s,t)= wW(j8;,s,f)|j8;|d|j8;|, (4.3.9) 
Jo 

where fij = \f}j\exp{i6). 

More details on the phase properties in terms of Pegg-Barnett formalism will be 
given in the following section. 



4.4 Phase properfies-Pegg-Barneff formalism 

Here we give essential background for Pegg-Barnett (55]- 1^ phase formalism. This 
formalism is based on introducing a finite (s + l)-dimensional space W spanned by the 
number states |0), |s). The physical variables (expectation values of Hermitian 

operators) are evaluated in the finite dimensional space W and at the final stage the 
limit s — > 00 is taken. A complete orthonormal basis of s + 1 states is defined on W as 



|e„) = -^=f^exp{ik@m)\k), {AAA) 

5 + 1 



where 



@m = ®0 + r. m = 0, 1 s. (4.4.2) 

s + 1 

The value of &o is arbitrary and defines a particular basis of s + 1 mutually orthogonal 
states. The Hermitian phase operator is defined as 

s 

^fre = 5Ii®m|e^)(e„|, (4.4.3) 

m=0 

where the subscript shows the dependence on the choice of ©o- The phase states 
l |4.4.1| are eigenstates of the phase operator l |4.4.3| with the eigenvalues @m restricted 
to lie within a phase window between @o and 27t + Oq. The expectation value of the 
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phase operator ( |4.4.3| in a pure state \xp) = Ylm=o ^m\j^)> where Cm is the weighting 
coefficient including the normalization constant, is given by 

s 

{M^e\^)= Y.@mm®m)f. (4.4.4) 

m=0 

The density of phase states is (s + 1)/(2;t), so the continuum phase distribution as s 
tends to infinity is 



P(0) = lim™^|(e^|^)|' 

(4.4.5) 

= ^E^V=oCn.C;,exp[f(m-m')e], 

where has been replaced by the continuous phase variable O. As soon as the 
phase distribution P{@) is known, all the quantum-mechanical phase moments can be 
obtained as a classical integral over @. One of the particular interesting quantities in 
the description of the phase is the phase variance determined by 

{{A^f) = j e2p(e)de - (^j ep(e)de^ . (4.4.6) 

As is well known the mean photon number and the phase are conjugate quantities in 
this approach and consequently they obey the following uncertainty relation 

((Acf.)2)((An)2)>l|([n,cf>])|2. (4.4.7) 

The number-phase commutator appearing on the right-hand side of l |4.4.7| can be 
calculated for any physical state |55]- as 



{[hA]) = i[i-27TP{@Q)]. (4.4.8) 

In relation to l |4.4.7| , we can give the notion of the number and phase squeezing |158l 
1159] through the relation 

S.^#^-l, (4.4.9) 
Sg = ^ - 1. (4.4.10) 
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The values of - 1 in these equations means maximum squeezing of the photon number 
or the phase. 

In this thesis we calculate the quantities given in this chapter for the models under 
discussion and then numerical simulations are performed using, e.g. Surfer or Matlab 
program. 
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Scientific results and their 
anal57sis: part I (static regime) 

In this part we investigate the quantum properties for a superposition of squeezed and 
displaced number states (SSDNS) without and with thermal noise. We suggest also a 
new type of multidimensional squeeze operator and discuss its properties in details. 

5.1 Quantum sfafisfical properties of superposition of squeezed 
and displaced number states 

This section is devoted to study properties of superposition of squeezed and displaced 
number states (SSDNS) such as orthogonality property, wave function, photon-number 
distribution, quasiprobability functions and phase properties in terms of Pegg-Barnett 
formalism. 

5.1.1 State formalism and some of its properties 

Here we develop a general class of quantum states, as a result of a superposition 
between two quantum states, described as a single mode vibration of electromagnetic 
field which is on a sudden squeezed-plus-displaced by a collection of two displacements 
180° out of phase given as 
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|r, a, n)e = Ae[£)(a) + e£)(-a)]5(r)|n), (5.1.1) 

where D(a) and S(r) are the displacement and squeeze operators, given by (|2.1.2|l and 
l |2.1.3|l . respectively, with ^1 = r and a are real; e = |e|e"^ is a parameter specified later 
on, (/) is a phase; |n) is the number (Fock) state, and Ae is a normalization constant 
given as 

|A2|~i = 1 + \ef + 2|e|exp(-2f2)Ln(4f2)cos(^, (5.1.2) 

where t = ae'', Ln(.) is the Laguerre polynomial. During the derivation of ( |5.1.2| , we 
used the relations 

^)^a)aD{a) = a + a, S'l'{r)dS{r) = dCj, - d^S^, ^ 
D(C)D(j8) = D(C + J8) exp[i(Cj8* - C*i8)], St{r)D{J3)S{r) = D(j8C, + j8*S,), 

and the relation [160] 

^exp(-f)P-^Lr^(f2) forn>fe, ^^^^^ 
f|exp(--|)(-f)*^""L^""(f2) fori?>n. 



{n\t){t)\k) = 



where = cosh r,Sr = sinhr; Ln(.) is associated Laguerre polynomial. The same 
states ( |5.1.1| were investigated independently in |161j . 

The density matrix p of this state can be written as 

P = |AeP(PM +P/)- (5.1.5) 

The part of the density matrix corresponding to the statistical mixture of two squeezed 
displaced number states is: 

Pm = D{a)S{r)\n){n\sHr)bHa) + \efb{-a)S{r)\n){n\sHr)bH-a), (5.1.6) 

while the quantum interference part has the form 

p, = e*b{a)S{r)\n){n\sHr)bH-a) + eb{-a)S{r)\n){n\sHr)bHa). (5.1.7) 

This quantum interference part of the density matrix contains information about the 
quantum interference between component states D(+a)S(r)|n) and this will be re- 
sponsible for some interesting behaviour of the phase distribution, as we will see. For 
completeness, the physical interpretation of such states can be related to a superposi- 
tion of coherent states formed due to two excitations on particularly excited harmonic 
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oscillators |[29]. It is clear that these states enable us to obtain generalizations of some 
results given in the literature, e.g. j70ll7 1 Il75ll88ll 1 62141 66] . Firstly, the wavefunction of 
state l |5.1.1t can be determined through ||3l] as 



Wn\x,r,a) = (x|r,a, n)e = 



9" 



[x\p.){ii\[D{a) + eb{-a)]S{r)exp{p.a^)\0) 



d'^li 



27T ' 

(5.1.8) 

where we have used the overcompleteness relation for coherent states; after substitut- 
ing for the element (x|^i) from |76] into j5.1.8| and using similar technique as in [31], 
we arrive at (with a is real) 



exp 



(5.1.9) 



+ eexp 



-^(x^+ y2a)2] Un[e'{x^^ + V2a)] 



where co and h are frequency of the harmonic oscillator and Planck's constant divided 
by 27r; they appeared here as a result of using the wavefunction (x |j8) when describing 
a coherent state of the harmonic oscillator in the coordinate representation in the 
derivation of l |5.1.9| . 

The inner product of the ket |r,a, n)^ with the bra e'(r',a',m| can be calculated 
using the wavefunction l |5.1.9|l through the relation 



e'(m,a',r'|n,a,r)e = 



dx'¥^m^*{x,r',a')¥n\x,r,a), (5.1.10) 

■oo 



together with the identity |167 



M /+^dxH^(x)H„(Ax + d)exp(-Mx2 + cx) = exp(^) 



M 



M-l\m/ M~A^\n\--iriin{m,n) jilml 



M I X]/=0 j\{n~j)\{m~j)l 



(5.1.11) 



2A 



t/(M-1)(M-A2) 



H 



H 



2^/{M~l]M J " ' \2^/M{M~^] J ' 



cA+2dM 



as 
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{m,a' ,r'\n,a,r) 



A^,Ae(IS!^)T ^min{m,n) nlm\[^^S^] 



VnlmlcoshR ^/=U ;!(n-;)!(m-;)! Vsinh 2R-' 

exp[|(tanhR - l)]H„_,-(Xi)H„_;(X2)[l + e'*e] 
+ exp[|(tanhR - l)]Hn_y(Vi)H^_;(V2)[e'* + (-l)("+")e] 



(5.1.12) 



and 



R 



r - r', Ti = (a - a')e^ T2 = (a + a')e^ Xi = 

, /<.;„i, on ' -^1 , /<.;„i, on ' 



(5.1.13) 

Vsinh 2R ' ""'^ ~ Vsinh 2fi ' ""^ ~ Vsinh 2R ' 

where Hjn(.) is the Hermite polynomial of order m. It is clear that for e' = -l,e = 1 
and n = m equation l |5.1.12t vanishes, that is the inner product of two states, one from 
even subspace and the other from odd subspace, is zero. Here we assume that e' = e. 
It is clear from l |5.1.12| that SSDNS are not orthogonal similarly as coherent states, but 
under a certain constrains they can be approximately orthogonal, e.g. when ti and 
T2 are very large such that tanhR < 1. It is easy to prove that for r' = r,a = a' and 
n = m equation ( |5.1.12| reduces to unity. 

When e = 0, a = 0, r = and r — > 0, l |5.1.12) reduces to the distribution coefficient 
for the superposition of displaced Fock states as 



(m|a, n, )e = AeVnlm! exp(-^) Y!,^o 



2^ •s^mln{m,n) (-l)""/g"+"'"2/ 
;!(n-;)!(m-;)! 



[1 + (-l)"+"'e] 



(5.1.14) 



A, 



-l)""Pa"+"'"2pexp(-^)[l + (-l)"+'"e]L^"''(a2), 



where p = min(n, m) and q = max(n, m). 

Further, taking e' = 0,a' = 0, and r' = in l |5.1.12|l and using the well-known 
definition of the photon-number distribution P(m) = \{m\r,a,n)^\'^, we get 



P'^Hm) = exp[f^(tanhr - 1)] Er^o 



•,min{n,m) n!m! 



Y r 2 1/r tanhr i^^ 



H 



(n~;)!(m -;)!;! 

2 



(5.1.15) 
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When r ^ 0, equation ( |5.1.15|l reduces to 




l + ]ep + ]e[ cos(j)e 



-a' 




mifi(fi,m) 




(5.1.16) 



X 



[(-!)""'■ + (-1) 



cos(/) + i|e|(-l)™ 'sin</)] 



It is clear that for e = and a = 0, P'^)(m) = 6n,m> as expected. Indeed, the os- 
cillations in the photon-number distribution of squeezed states [168| are important 
effect which emphasizes the deviation from the Poissonian distribution. These oscilla- 
tions are remarkable for large squeezing parameter and can be interpreted using the 
Bohr-Sommerfeld band picture in phase space. Here we study the photon-number 
distribution for SSDNS with particular attention to the effect of the overlap between 
oscillators. During our numerical analysis we use only three values for e, namely, 1,-1 
and i, respectively in correspondence to even squeezed and displaced number states 
(ESDNS), odd squeezed and displaced number states (OSDNS) and Yurke squeezed 
and displaced number states (YSDNS). These states have been considered in connec- 
tion with, for r = 0, n = 0, even coherent, odd coherent and Vurke-Stoler states. 



Figure 5.1: Photon-number distribution P(m) of ESDNS (solid curve) and YSDNS (dashed curve) for 
a = 2, n = 5 and r = 0.8. 

From Fig. 5.1 it is clear the oscillatory behaviour of squeezed states as expected. 
Comparison of the solid and dashed curves shows that interference part (the overlap 
between different oscillators, i.e. p/) increases the oscillations in the photon-number 
distribution. Intuitively, when n increases the oscillations in P(m) become more pro- 
nounced as the zeros of the Hermite polynomial increase. Furthermore, this situation 




m 
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is still valid with increasing r f68]. It is convenient to point out that the numerical 
analysis of P(m) for OSDNS has a similar behaviour as for ESDNS, whereas for the 
statistical mixture it is similar to that of YSDNS, i.e. it does not depend on the quan- 
tum interference between the components of the states D{+a)S{r)\n). We conclude 
this subsection by discussing the sub-Poissonian behaviours of SSDNS with the help 
of g'^HO). This can be done by calculating the expectation values of a^a and d^d^ in 
terms of the state l |5.1.H and inserting the results into ( |4.1.1| (with t = 0). The cal- 
culations of the previous moments are rather lengthy but straightforward, where the 
relations ( |5.1.3| and l |5.1.4| should be frequently used. 

a) b) 
C C 




Alpha Alpha 
(a) (b) 

Figure 5.2: Normalized normal second-order correlation function gf'^'(O) a) for VSDNS, n = and r = 0.5 
(solid curve), 0.9 (dashed curve), the straight line is corresponding to g'^' (0) for the coherent state; b) for 
VSDNS, n = 4 and r = (solid curve), 0.5 (long-dashed curve), 0.9 (short-dashed curve), the straight line 
is corresponding to g'^' (0) for the coherent state. 

Here, we do not write down their expressions. Now we start our numerical ana- 
lysis with particular attention paid to Vurke squeezed and displaced number states. In 
fact, Vurke and Stoler |165] have shown that a coherent state propagating through 
an amplitude dispersive medium, under specific conditions of parameters, can evolve 
into a superposition of two coherent states 180° out of phase. For such Yurke-Stoler 
states g'^HO) = ^> which corresponds to Poissonian distribution. As soon as r in- 
creases (r = 0.5,0.9), the Poissonian distribution disappears and g'^^O) takes on values 
corresponding to super-Poissonian statistics in the short starting interval of a and 
sub-Poissonian behaviour elsewhere (Fig. 5.2a). Also it is clear that the length of the 
sub-Poissonian interval decreases as the value of r increases. On the other hand. 
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when the number of quanta increases (e.g. n = 4) for the same values of r, we see 
that g'^^O) takes on always super-Poissonian values for large value of r (r = 0.5,0.9) 
and it has sub-Poissonian behaviour when r = in the short starting interval of a (Fig. 
5.2b). 

In general we noted g^^\0) for ESDNS and OSDNS has the same behaviour as 
for YSDNS. In other words, the known behaviours of g'^^O) for even coherent states 
and odd coherent states are changed by increasing the values of r and n, and the 
sub-Poissonian interval decreases as r increases. 

5.1.2 Quasiprobabilify distribution functions 

Here we examine the quasiprobability distribution functions for SSDNS. We can adopt 
two types of these functions which are Q-function (Husimi) and W-function (Wigner). 
We use the s-parametrized characteristic function C(C, s) as defined in section 4.3 to 
calculate such functions. Using the definition of p given by l |5.1.5| together with the 
relations l |5.1.3| -( |5.1.4| in a straightforward way, we get 

C(^)(Cs) = |A,|2exp(-|fe2|/2 + ||^|2)L„(|ij2|) [e^p[(^ _ 



+ exp[(C - CMef] + |e|A.|2 [exp(-i(/, - ^ + ^\Cf)U{\kl\) 



(5.1.17) 



where 



+ exp(i(/,-^ + ||C|2)L„(|^? 



k = CCr + C*Sr, k+ = {C + 2a)Cr + (C + 2a)S^, 
J?_ = (C - 2a)Cr + iC - 2a)Sr. 



(5.1.18) 



As soon as the characteristic function is obtained, the W-function and Q-function 
can be calculated easily by inserting j5.1.17| into j4.3.4| and carrying out the integration 
(see Appendix A for the details), we get 



W(e)(x,y) = ^tllM. {|e|2exp[-2(y2e-2'' + e^''{x + a)2)]Ln[4(y2e-2'' + e^''{x + af)] 

+ exp[-2(y2e"2'' + e^''{x - af)]Ln[A{y^e-^'' + e^''{x - af)] (5.1.19) 

+2|e| exp[-2(y^e"2'' + e2''x^)]Ln[4(y2e"2'' + e^'^x^)] cos(4ya -(/))}, 
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|Aj[2(!S!^)" 
7tn \ coshr 



{exp[-|]8 - af- - Re(]8 - a)^tanhr] 



^ J^]!' + exp[-|(|i8 - af + |j8 + af)] 
Re [exp[-^si*i:[(j8 - a)2 + (j8 + a)^] + 2ialmj8 + 



(5.1.20) 



i(j3''+a) 



]Hn[: 



- Vsinh 2r Vsinh2r-' 

where ]8 in j5.1.19| has been taken as x + iy and x and y which are corresponding to 
the quadratures. 

The interference in phase space which is representative to the superposition states 
may be seen clearly in the behaviour of the W-function. In general we find that 
there are two regimes controlling the behaviour of W-function as well as the phase 
distribution for the states (|5.1.1t which are a >> 1 and a < 1 provided that n is 
finite. For the first regime, the even- and odd-cases (i.e. e = 1 and -1) provide 
similar behaviours and there is one-to-one correspondence between the components 
of density matrix l |5.1.5| and that corresponding to the curves. Nevertheless, all these 
features are washed out in the second regime and the behaviour dramatically change. 




(b) 



Figure 5.3: W-function against r and y for n = l,a = I and for a) e = 1; b) e = ~1. 



More illustratively, W-function j5.1.19| includes three terms, the first two terms rep- 
resenting the statistical mixture of squeezed displaced number states and the third one 
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is the interference part. When a is small (as shown in Fig. 5.3a and b for superposition 
of displaced number states) the contributions of these components are comparable so 
that even- and odd-cases distributions have different shapes. On the other hand, when 
a is large the structures of the W-function for even- and odd-cases are almost similar, 
i.e. they include two symmetrical peaks originated at +a and interference fringes are 
in between. Of course, such situation is still valid if squeezing in the displaced super- 
imposed number state in optical cavity is considered, however, the peaks will be then 
stretched. This behaviour is reflected in the behaviour of the phase distribution, as we 
will see. 

The quantum-interference term is not visible in the Q-function jl69j . For ESDNS 
the Q-function is shown in Figs. 5.4a-d for shown values of the parameters. By 
numerical analysis we obtained for r = and n = two identical Gaussian shapes 
which are stretched by increasing the value of r. As n increases (n = 1), having r = 0, 
the peaks for r = 0, n = are replaced by two identical closed top hole peaks, each of 
them represents the Q-function of the displaced Fock state and there is evidence 
for splitting in the internal edges of the holes in Fig. 5.4a. Increasing r(r = 0.9) for 
n = 1 as before, the two peaks for r = are completely separated from each other, by 
the action of squeezing, and each of them is split to two identical Gaussian-like shapes 
with the same root (Fig. 5.4b). However, when n = 4 and r = 0.5,0.9, we observe 
that for r = 0.5 the two peaks occurred for n = 1 are replaced by two symmetrical 
stretched circular holes (Fig. 5.4c). But for r = 0.9 each hole is converted to two 
symmetric peaks separated by four spikes (Fig. 5.4d). 

The Q-function for OSDNS and YSDNS can be represented by similar figures. This 
shows that the interference term plays negligible role in the Q-function. 



5.1.3 Phase properties 

We shall use the relations of the Pegg-Barnett formalism given in section 4.4 to study 
the phase distribution for the superposition of displaced and squeezed number states 
( |5.1.1| . In this case Cm in l |4.4.5| is given by Cm = Cm(r,a,n, e) = (m|r,a,n)e and this 
quantity can be deduced from l |5.1.12| by suitable choice of the parameters. 

The photon-number variance is defined as 
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Figure 5.4: Q-function for ESDNS and a = 2 for a) n = 1 and r = 0; b) n = 1 and r = 0.9; c) n = 4 and 
r = 0.5; d) n = 4 and r = 0.9. 



((An)2) = (n2) - (n)2 

(5.1.21) 

= (d7^d2) + (n) -(n)2, 

where the number operator n = d^d. In terms of l |5.1.1t the quantities (d^d^) and (n) 
can be straightforwardly calculated. For completeness, the phase variance of l |5.1.H 
reads 

((Aci.)2) = + 4ReE„,>n,'C„(r,a,n,e)C;;,(r,a,n,e)tl);^ 

(5.1.22) 

I^ei Em>m' Cm{r,a,n, e) C;;, (r , a, n, e) 
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The value :jt^/3 is the phase variance for a state with uniformly distributed phase, e.g. 
the vacuum state. 

At the mean time, for understanding the behaviour of the phase distribution of the 
states l |5.1.1t it is more convenient firstly to study such properties for two subsidiary 
states which are the superposition of displaced number states fl70j, and squeezed and 
displaced number states |29J. These two states are of interest because it has been 
shown for the former states that they can exhibit strong sub-Poissonian character as 
well as quadrature squeezing. The bunching and antibunching properties for the latter 
states have been discussed in |30]. Investigation of such properties for these two kinds 
of states most clearly illustrates the role of the different parameters in the state jS.l.H 
in the behaviour of the distribution. In the following we discuss the phase distribution 
determined by ( |4.4.5| , and the phase variance and phase squeezing, respectively. For 
simplicity we restrict our investigation to real values of a and e. 

(a) Phase probability distribution 

As we mentioned before, there are two regimes controlling the behaviour of the phase 
distribution for the states under discussion similar to the W-function. We start our 
discussion by focusing the attention on the behaviour of the superposition of displaced 
number state. It is important to mention that the properties of displaced number states 
D{a)\n) have been given in |166) and their phase properties in [ITlj. In fact the phase 
properties of such states are interesting since they connect the number state, which 
has no phase information, and coherent states, which play the boundary role between 
the classical and nonclassical states and always exhibit single peak structure of phase 
distribution. The phase distribution of displaced number states exhibits nonclassical 
oscillations with the number of peaks which is equal to n + 1. This result is interpreted 
in terms of the area of overlap in phase space jl71] . Unfortunately, the situation is 
completely different for the superposition of displaced number state (see Figs. 5.5a 
and b for shown values of the parameters). From Fig. 5.5a it is clear that the quantum 
interference between component states D{a)\n) and D{-a)\n) leads to the lateral non- 
classical oscillations which become more pronounced and narrower as n increases. 
However, the statistical mixture of displaced number states shows n + 1 peaks around 
= 0. We have not presented the phase distribution for the odd-case since it is sim- 
ilar to that of the even-case. Comparing this behaviour with that of even (or odd) 
coherent states, we conclude that the number states in the superposition make the 
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phase information more significant (46J. Turning our attention to the Fig. 5.5b where 
a = 1, we can see that the distribution of displaced number states exhibits two-peak 
structure as expected for n = 1 |171| , however, the distribution of even-case displays 
one central peak at = and two wings as ^ +7t, and finally the distribution 
of odd-case provides four peaks. That is the distribution here is irregular, however, 
more smoothed than before and the structure of the statistical mixture (central part) 
is modified by the action of the interference term arising from the superposition of 
the states; this was the case for W-function. 



b) 



a) 




Theta 
(a) 




Theta 



(b) 



Figure 5.5: Phase distribution P{&) for superposition of displaced number states and for a) a = 2, e = 1 
, n = 1 (dashed curve), and 2 (solid curve); b) a = 1, n = 1 and e = (bell-centered curve), 1 (solid curve), 
and -1 (dashed curve). 



Before discussing the phase properties of the displaced and squeezed number states 
it is reasonable to remind the behaviour of the well known squeezed states. As known 
for squeezed states with non-zero displacement coherent amplitude, the phase dis- 
tribution exhibits the bifurcation phenomenon. In this phenomenon the single peak 
structure of the coherent component is evolved into two peaks structure with respect 
to both a (for large fixed value of squeezing parameter r) and r (when a takes fixed 
value) 160]. This phenomenon has been recognized as a result of the competition 
between the two peaks structure of the squeezed vacuum state and the single-peak 
structure of the coherent state. For squeezed and displaced number states such phe- 
nomenon cannot occur due to the effect of the Fock state which replaces the initial 
peak (r = 0) for coherent state by a multi-peak structure, i.e. by n + 1 peaks (see 
Fig. 5.6a for shown values of the parameters). From this figure one can observe that 
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there is a three-peak structure corresponding to the case n = 2 of displaced Fock 
state. The height of the central peak (i.e. at = 0) is almost the same and equals 
{^/^^)\Yl'm=o^mf- That is the central value of the phase distribution P{@ = 0) is in- 
sensitive to squeezing provided that r is finite. However, the lateral peaks undergo 
phase squeezing as r increases, i.e. the peaks become narrower. 

Now we can investigate the behaviour of the superposition of displaced and squeezed 
number states (see Figs. 5.6b-d for shown parameters). Figs. 5.6b and c are given for 
the second regime for even- and odd-cases, respectively. From Fig. 5.6b we can see 
that the initial oscillations are increased compared with Fig. 5.6a as a result of the 
interference in phase space. Further, the initial lateral peaks can evolve in the course 
of increasing r to provide bifurcation shape, i.e. the distribution curve undergoes a 
transition from single- to a double-peaked form with increasing r; however, the central 
peak is almostly unchanged with increasing r. This peak splitting is connected with the 
squeezed states ||60j. Further, as the number of peaks increases for r > 0, the distribu- 
tion becomes more and more narrower. On the other hand, the phase distribution of 
the odd-case is quite different as we have shown before where the initial peaks are not 
significantly changing as in the even-case. More precisely, initial distribution becomes 
broader for a while and suddenly (at r ~ 0.5) breaks off to start to be a narrower 
distribution for later r. It is clear that in this regime the state l |5.1.1| becomes more 
and more nonclassical and the even- and odd-cases are distinguishable. Further, the 
phase distribution of the even-case is more sensitive with respect to squeezing than 
the odd-case. Nevertheless, for the first regime where a is large we noted that the 
phase distribution carries at least the same initial information regardless of the value 
of r (see Fig. 5.6d for the shown values of the parameters). 
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Figure 5.6; Phase distribution P{&) for superposition of displaced and squeezed number states for 
{e,a,n) = (0, 1,2), (1, 1,2), (-1, 1,2) and (1,2,2) corresponding to the cases a,b,c and d, respectively. 
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(b) Variance, amplitude and phase squeezing 

Here we investigate the behaviour of the phase variance, and amplitude and phase 
fluctuations following l |5.1.22| , ( |4.4.9| and j4.4.10| , respectively. We start our discussion 
by analyzing the behaviour of the superposition of displaced number states. For this 
purpose Figs. 5.7a and 7b are shown for the phase variance, and amplitude and phase 
fluctuations, respectively. It is seen that in general the phase variance starts from the 
value :jt^/3 (the vacuum state value) and returns back to it when a is large, but through 
different routes. To be more specific, the phase variance of displaced number states 
starts from the value for the vacuum, goes to a minimum, and then comes again to 
7r^/'5. However, the behaviour of the superposition of displaced number states takes 
different ways to arrive at the same result, i.e. it starts from ir^/'i as before, goes to 
the maximum value and eventually comes back to the value of vacuum. The compar- 
ison of the two cases shows the role of the quantum mechanical interference between 
state components. Further, as n increases, the oscillations in the variance become 
more pronounced. Comparison of the behaviour of even- and odd-cases shows that 
they are different only over the initial short interval of a, i.e. when a is small, and 
this agrees with what we have discussed earlier. So we can conclude that for intens- 
ities high enough of the coherent field, the variance of the phase is approximately 
randomized. The route to this randomization is dependent on the choice of e. With 
respect to the amplitude and phase fluctuations, we can note from j4.4.8| that these 
quantities depend not only on the intensity of the field, but also on the choice of the 
reference angle @o. We have chosen here @o = -tt, where the mean value |([iV,4>0])| 
of squeezed displaced number state approaches unity. Fig. 5.7b has been obtained to 
illustrate the parameters Sm and Sg which provide information about the degree of 
squeezing in n and ^g. One can observe from this figure that when e = (displaced 
number state) and a — > 0, the parameter Sjv tends to - 1, which means that the number 
state is 100% squeezed with respect to the operator n. This situation is expected since 
((An)^) = for the number state. Further, the larger the number of quanta is the 
shorter is the interval over which Sm is squeezed. Also when a >> 1, squeezing in 
Sg is remarkable, whereas Sm becomes unsqueezed. This result can be deduced from 
the behaviour of the phase variance (see Fig. 5.7a), where ((A<I>)^) ~ at a 2^ 4 and 
this should be connected with maximum squeezing in Sg at this point. Such behaviour 
of Sm and Sg confirms the fact that the number of photons and phase are conjugate 
quantities in this approach. As is known displaced number states are not minimum 
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uncertainty states and the variances for the quadrature operators never go below the 
standard quantum limit. Moreover, they may exhibit sub-Poissonian statistics for the 
range < 1/2 |166| . However, there is no relation between the sub-Poissonian stat- 
istics and the fluctuation in the amplitude or the phase. This fact has been shown 
before for the down-conversion process with quantum pump where the signal mode 
can exhibit amplitude squeezing and at the same time it is super-Poissonian |172| . 

We proceed by discussing the behaviour of the superposition states (long-dashed 
and circle-centered curves) in Fig. 5.7b where the interference in phase space starts 
to play a role. We noted (from our numerical analysis) that only the even case can 
provide squeezing in Sm with maximum value at the origin and squeezing interval 
larger than that discussed before. Indeed, this maximum value is related also to that 
of number state at a = 0. It should be stressed here that a = for the odd-case may 
lead to a singularity. As we mentioned before the superposition of displaced number 
states can exhibit strong sub-Poissonian character as well as quadrature squeezing 
|170| . Now we can illustrate the role of the squeezing in the superimposed displaced 
number states optical cavity with respect to variance, amplitude and phase fluctuations. 
It is obvious, when squeezing is considered, that the initial value (at a = 0) of the phase 
variance is shifted since we have initially squeezed number state which is providing 
phase information. However, when a is large and r is finite or also r is large, it 
can be proved simply that the coefficient Ciniria.nie) vanishes and consequently the 
phase variance tends to ;t^/3 (becomes randomized). Moreover, the routs are here 
similar to those of Fig. 5.7a. On the other hand, squeezing could be seen only in Sg 
for e = (see squared-centered curve in Fig. 5.7b). Furthermore, comparison of the 
short-bell-centered curve (of displaced Fock state) and squared-centered curve reveals 
that squeezing parameter reduces the amount of squeezing in Se, too. This means 
that the superposition of displaced and squeezed number states provides quadratures 
squeezing which possesses less information about amplitude and phase fluctuations. 

5.1.4 Scheme of production 

Nonclassical states of light in cavity are of increasing importance in quantum optics. 
There are two principal approaches how to generate them [35lll73j : (i) Find the 
appropriate Hamiltonian which transforms via unitary time evolution the initial states 
to the desired final state, e.g. Schrodinger cat states, (ii) Make a measurement on 
one of two entangled quantum systems and obtain the state of the other system by the 
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Figure 5.7: a) Phase variance of a superposition of displaced number states against a for (e,n) = (0,3) 
(long-dashed curve), (1,1) (bell-centered curve), (1,3) (solid curve) and (-1,3) (short-dashed curve). The 
solid straight line is corresponding to the phase variance of vacuum, b) Amplitude and phase fluctuations 
of a superposition of displaced and squeezed number states against a for (e, n, r) = (0, 1, 0) ( solid curve 
and Sg short-bell-centered curve), (0, 2, 0) ( Sjv short-dashed curve and Sg long-bell-centered curve), (1,1,0) 
(Sn long-dashed curve), (1,2,0) (Sjv circle-centered curve) and (1,2,0.5) (Se squared-centered curve). The 
solid straight line is the bound of squeezing. 



corresponding state reduction. 

We use here, to suggest generation of superposition of such states, quantum state 
engineering, which is based on injecting JV two-level excited identical atoms propagat- 
ing through vacuum cavity field. The atoms passing the cavity one by one in such a 
way that each injected atom increases the number of Fock states, are building up the 
cavity field state by one. In microwave experiments this is achieved by using Rydberg 
atoms and very high-Q superconducting cavities so that spontaneous emission and cav- 
ity damping are quite negligible on the time scale of the atom-field interaction. The feth 
atom enters the cavity in a given coherent superposition state Cg\g)k + C^\e)k, where 
\g) (|e)) denotes ground (excited) atom state and Cg and Cg are the corresponding su- 
perposition weights, and it can be detected in the ground state. The atom-field resonant 
interaction is described by Jaynes-Cummings model via the interaction Hamiltonian 
H = hx{do+ + d^cr_), where cr+ and cr„ are atom raising and lowering operators and 
X is atom-field coupling constant. Once the JV atoms crossed the cavity, leaving their 
photon in the cavity and decaying to the ground state, the pure field state in the cavity 
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is Hzi] 



N 



(5.1.23) 



m=0 



where A is a normalization constant and the coefficients come from the old coef- 
ficients A^~^ and A'^^J. via the recurrence formula 



= (1 - SmMllWe smixT^VJ^) + (1 - 6rn,N)At'ci^^ cos{xTnVJ^), (5.1.24) 



where Aq = 1 corresponds to the initial vacuum state and Tjv is the interaction time 
of the JVth atom with the field. 

Now SSDNS can be represented in the Fock state basis as 



where Cm = (m|r,a,n)e is the expansion coefficient, transition amplitude, and can 
be obtained from l |5.1.12|l . Then by controlling Cg, Cg, the atomic superposition coef- 
ficients, we can prepare each atom, in principle, before entering the cavity in such a 
way that A^ = Cm- It must be held JV > n, where n is the mean number of photons 
for the desired state. In this way the desired state may be obtained. The solution of 
the recurrence relations j5.1.24| , that give the unknown coefficients Cg, Cg in terms of 
the known coefficients Cm, is available [5511175] . 

Recently, a great progress has been done in quantum states generation using trap- 
ping ions and laser cooling (55|. An ion trapped in a harmonic potential can be 
regarded as a particle with a quantized center-of-mass motion. One can consider an 
ion trapped in a harmonic potential and driven by two laser beams tuned to the first 
lower and upper vibrational sidebands, respectively jl75] . The author showed that, in 
the rotating-wave approximation, considering sideband limit, i.e. the vibrational fre- 
quency is much larger than the other characteristic frequency of the problem, and the 
behaviour of the ion in the Lamb-Dicke regime, i.e. 77 « 1, with choosing the amp- 
litudes and phases of the lasers in the appropriate forms, the interaction Hamiltonian 
describing the problem takes the form 



where g = riQexp{-^), and S_,S+ are the raising, lowering operators for the two- 
level ion. Further, Q is the common Rabi frequency of the respective lasers and 77 is 




(5.1.25) 



m 



Hi = -ig(d + d^)(S_ + S+), 



(5.1.26) 
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the Lamb-Dicke parameter which is connected with the vibrational frequency v, the 
wave vector of the driven field k and the mass of the trapped ion M are related by 
77 = Assuming the ion is initially in the ground state \g) and the motion is in the 

arbitrary state \^), then after an interaction time t the system is in the entangled state 

|i/)(t)) = cos(gT6)|(/))|g) - isin(gTO)|(/))|e), (5.1.27) 

where O = -i(d - d^). Detecting the internal state of the ion and considering the ion 
in the ground state, then the vibrational motion collapses to 

\^)g = Ng[b{a) + D(-a)]|(/)), (5.1.28) 

where a has been considered equal to gr, D{a) is the displacement operator as before 
and JVg is a normalization factor. If \^) is a squeezed number state f7l ll65j we get the 
motional even squeezed and displaced number state of ion. On the other hand, if we 
find the ion is in the excited state similarly we can get odd squeezed and displaced 
number state. 

Furthermore, the quantum states discussed in this section can be generated from 
a micromaser, using quantum superposition of squeezed and displaced states from 
nonlinear optical processes, such as optical parametric generation and amplification, 
to which a radiation in the Fock state is initially introduced. 

5.1.5 Conclusions 

Superposition of quantum states controls the behaviour of single states by making it 
more or less pronounced, as a consequence of the interference between different prob- 
ability amplitudes. We have introduced a superposition of two squeezed displaced num- 
ber states and discussed the quantum statistical properties of the resulting state. We 
concentrated on three special sub-states, ESDNS, OSDNS and YSDNS corresponding to 
e = 1, -1 and i, respectively. For these states we studied scalar product, photon-number 
distribution P(m), normalized second-order correlation function g'^^O), W-function, Q- 
function, and phase properties. By means of the scalar product of SSDNS we have 
demonstrated that the states which belong to the same subspace are not orthogonal 
in general, but if r\ and T2 are large such that tanh f? < 1, they become approximately 
orthogonal. The oscillations in P(m) have been demonstrated. These oscillations are 
more pronounced as r and n increases. Such oscillations have been recognized as 
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a striking feature of highly nonclassical states. The normalized intensity correlation 
function g^'^\0) for this superposition shows that the sub-Poissonian interval decreases 
as r increases. For quasiprobability distribution functions for n f but finite, we have 
noted that for W-function there are two regimes controlling its behaviour depending 
on whether the superposition is macroscopic (a >> 1) or microscopic (a < 1). In the 
first regime, the even- and odd-cases (i.e. e = 1 and -1) give similar behaviours and 
the structure of the density matrix is remarkable in figures. In other words, for this 
case Wigner function includes two symmetrical peaks originated at +a and interfer- 
ence fringes are in between. However, for the second regime the structures of the 
W-function for even- and odd-cases are quite different. Further, in general we noted 
that there are negative values of the W-function for different e, which are signature 
of their quantum mechanical features. For all cases squeeze parameter is respons- 
ible for peaks stretching, which increases as r increases, and the occupation number 
n is responsible for "chaotic" behaviour, which is more pronounced as n increases. 
The Q-functions exhibit a number of interesting quantum effects. For example, the 
Q-function of shifted ground state, e.g. a f 0, r = 0, n = 0, is squeezed by increasing 
the value of r and has several peaks and spikes by increasing the number of quanta n. 
Also it is noted that the quantum-interference term is more visible in the W-function 
but not in the Q-function [169]. For the phase distribution we noted that in the first 
regime, the even- and odd-cases give similar behaviours and the structure of the dens- 
ity matrix is remarkable in figures similar to W-function. All these facts are washed 
out in the second regime where the behaviour becomes irregular, however, smooth. 
In general we noted that the higher the number of quanta is, the more peaks the 
distribution possesses. Influence of the squeezing parameter could be recognized in 
the second regime where the distribution exhibits peak-splitting and peak-narrowing 
and this is in contrast with the fist regime. For the phase variance we conclude that it 
asymptotically goes to the value 7t^/'5 of the uniform distribution when either a or r 
is large, but through different routes. We have shown also that this superposition can 
exhibit photon-number fluctuations and phase fluctuations. 

Finally, based on recent results [17411175] we have suggested ways for generation 
of such superposition in the framework of micromaser [174] and trapped ions |175] . 
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5.2 Quantum sf at isf ical properties of superposition of squeezed 
and displaced states with thermal noise 

In this section we investigate the influence of thermal noise on the properties of the 
previous states l |5.1.H using density operator formalism. 

5.2.1 Density operator for the states 

Sometimes we have not enough information to specify completely the state of the 
system and hence we cannot form its wave function and consequently the density 
operator must be used to describe the state of the quantum-mechanical system. The 
well-known example for such a state is provided by the thermal excitation of photons 
in a cavity mode maintained at the temperature T, which is described by the density 
operator, pch, as |176| 



where z = is the quotient of Bose-Einstein (geometric) distribution and mean 

number of photons n = ^exp(jl^) - 1 j for thermal equilibrium at temperature T, v 
being radiation frequency, Kb is the Boltzmann constant and h has the same meaning 
as before. It is clear that the thermal distribution has a diagonal expansion in terms of 
the Fock state |n). This diagonality causes the electric field expectation value to vanish 
in thermal equilibrium. Such thermal field can be generated by a thermal source 
composed of many independent atomic radiators and consists of the superposition of 
waves of many different frequencies within some continuous range. These waves can 
be regarded as independent waves with random phases [88]. 

In the following we define the thermal superposition squeezed and displaced states 
(TSDS) similar to that of squeezed thermal state |70]- 173], displaced squeezed thermal 
state |74] and a superposition of coherent state with thermal noise jl77] , as a sum of 
states jS.l.H weighted by the Bose-Einstein distribution, which can be represented by 
the density operator in the form 




(5.2.1) 



oo 



Pt = 



1 + n 



^z"[D(a) -h eb{-a)]S{r)\n){n\S{-r)[b{-a) + e*D(a)], 



(5.2.2) 
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where all the notations have the same meaning as before and A^- is the normalization 
constant, obtained with the aid of the normalization relation Trp = 1 to take the form 



It is noted from equation ( |5.2.2| that the thermal superposition of displaced and 
squeezed number states which has been considered in optical cavity, is not a product 
of the states for independent oscillators. In fact we can find interesting nonclassical 
effects which can occur as a result of existing correlations between oscillators; this will 
be seen later. 

As before the parameter e can take either the value 1, or -1 or i corresponding to 
thermal squeezed even displaced state (TSEDS), thermal squeezed odd displaced state 
(TSODS) and thermal squeezed Yurke displaced state (TSYDS), respectively. 

In the following we calculate the quasiprobability distribution functions which will 
be used to derive the second-order correlation function, photon-number distribution, 
purity, and phase distribution. 

5.2.2 Quasiprobabilify distribution functions 

Using the s-parametrized characteristic function given in section 4.3 together with 
(|5.2.2t we straightforwardly get 



X^^ = 1 + \ef + 2|e| cos (j) exp[-2(l -h 2n)t^] 



-1 



(5.2.3) 



C(C,s) = A^^^exp [-(A -h i[l - s])|C|2 - f (C^ + C')] 
X |exp(Ca* - C*a) + |e|2exp(C*a - Ca*) 

+ |e|exp[-i(/) - G - (CC* + C*C)] + \e\exp[i(}) -G + (CC* -h C*C)] 



(5.2.4) 



where 



G = (1 + 2n)[2|a|2cosh2r + (a^ + a*2)sinh2r], 
C = (1 + 2n) [a cosh 2r + a* sinh 2r] , 
A = n + (1 + 2n)S2, E = (1 -h 2n)SrCr. 



(5.2.5) 
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During the derivation of l |5.2.4| , relations ( |5.1.4| and (A.2) in appendix A have been 
used. Consequently, using l |4.3.4| together with the identity (A.4) in appendix A, the 
s-parametrized quasiprobability distribution functions are 



W(j8,s) = 



X J. exp 



+ |epexp 



(A+i[l-s])2-E2 



[A+^{i~s)]\p+af + ^E[{p+af + {p*+a*f] 



(5.2.6) 



(A+^[l-s])2-E2 



+ |e| [exp - G + K) + exp {-i(p - G + K*)] , 



where K is given by 



_ [A + i(l - s)](C - P){C* + p*) - E[(ReC - ilmpf + (ilmC - Re^f] 

(A + i[l - s])2 - E2 ' I • • J 

where Re and Im mean real and imginary parts. As we mentioned before when 
s = 1,0, and -1, expression j5.2.6| gives Glauber P-function, Wigner W-function and 
Husimi Q-function, respectively. 

By taking the appropriate values of the parameters e, a, r, one can easily use ex- 
pression l |5.2.6| to obtain the corresponding ones for some special states, such as 
thermal squeezed coherent state [74], and squeezed thermal state 170] [71], also for 
thermal odd- and even-coherent states |177| . 

In the following we examine the Q-function (s = -1), and the W-function (s = 0) 
for TSVDS given by l |5.2.6t . It should be noted that for all the cases we assume the 
parameter a is real and equal to 2. In Fig. 5.8 we have plotted the Q-function against 
both X = Rej8 and y = Imj8, for r = 0.5, and n = 2.5. From this figure we can observe 
two identical peaks, both of them are stretched under the action of the squeezing 
parameter r. Further, these peaks are broader than those of the statistical mixture of 
coherent states and this is connected with mean number of thermal photons, i.e. for 
large n the peaks become more and more broader. 
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Figure 5.8: Q-function for TSYDS and a = 2 and r = 0.5 for n = 2.5. 
a) b) 




(a) (b) 
Figure 5.9; W-function for TSVDS and a = 2 and r = 0.5 for a) n = 0.5 ; b) n = 1. 

In Figs. 5.9a and b we have plotted the W-function against x = Rej8 and y = Imj8 
for r = 0.5, n = 0.5 and for n = 1, respectively. In Fig. 5.9a we can see the W-function 
in X - y plane, showing that the oscillator phase space consists of two Gaussian bells, 
corresponding to statistical mixture of individual composite states; also we can see 
interference fringes inbetween arising as a result of the contribution originating from 
the quantum superposition. Furthermore, we find that the W-function takes some 
negative values, which is indicating nonclassical effect, due to the interference of the 
states in phase space. Increasing thermal photon number n and keeping the value of 
r as before, we can realize that the two Gaussian bells are stretched in the y-direction 
and the interference fringes become more sharper than before. In the mean time the 
negative values in W-function start to smooth out reflecting the classical effect of the 
thermal noise in the system (see Fig. 5.9b). 
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5.2.3 Second-order correlation function 

Here we study the sub-Poissonian statistics for the state under consideration. So that 
the relation \L'5.5\ has been used to deduce the expectation values (d^d) and (d^d^) 
as 



(d^d) = A^^' 



(A + |a|2)(l + \ef) + 2|e|e"^(A - \Cf) cos(j) 



(5.2.8) 



+2|e|e-^ [2A2 + - A\CfA + |C|^ - E{C^ + C*^)] coscp}. 



(5.2.9) 




Figure 3.10: Normalized normal second-order correlation function gr'^'(O) for TSODS for n = 0.1. 



As is well-known in the thermal optical cavity, photons have tendency to bunch 
each other, where thermal light can be described as a classical light. In the thermal 
coherent squeezed optical field in cavity the loss of nonclassical effects is much faster 
than at zero temperature. This is noticed from the numerical investigations of the 
function g^'^^{0), for both TSEDS and TSVDS; for example when we take n f 0, we find 
the values of g'^^O) always greater than one showing super-Poissonian behaviour 
whatever the values of r and a are. On the other hand, for TSODS, and for small 
value of n we realized that there is sub-Poissonian behaviour with maximum value 
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when a and r are small, which indicates that the states are closed to the odd coherent 
states, see Fig. 5.10. Finally we would like to point out that for n > 1 sub-Poissonian 
behaviour entirely disappeared. 

5.2.4 Photon-number distribution 

Here we discuss the photon-number distribution P(n) for the states given by l |5.2.2| . 
As we mentioned before, this function can be deduced using the relation between 
Laguerre polynomial and Wigner function given by j4.3.8| . In this case the integration 
may be carried out applying the same technique as in Appendix A and the photon- 
number distribution reads 



Pin) = 



V(A+1)2-E2 



(1 -h |e|^)L„^(-X_)L;f„(-X+)exp 



[(A+l)[a|'^ + E(af-a|)] 



(5.2.10) 



+2|e|L^^(C2v„)L;L(C2y+)exp 



[{A+i)\Cf~E{Cf~Cl)] _ ^ 
(A+l)^^-E2 *J 



COS (f) \ , 



where 



v+ = [A + l± E]~K X+ = al[{A + E){A + E+ 

X_ = af[(A - E)(A - E + 1)]"^ V± = [(A + E)(A + E + 1)]-^, 

and Ln(.) is again the associated Laguerre polynomial. 

Now when we take r — > in (|5.2.10t and using the identity |!178] 



(5.2.11) 



(5.2.12) 



we arrive at 



Pin) ^ i^r 



+2|e| exp 



(1 + |e|^)exp 



■ l+n J I l+n 



(2n + l)|a|2 
l+n 



(5.2.13) 
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which is in agreement with equation (26) of |177] provided that we set |e| = 1 for even- 
and odd-thermal coherent states. 



0. 10q 




0.00 1 

0.00 5.00 10.00 15.00 20. 00 25. 00 

n 

Figure 5.11: P(n) for TSEDS for r = 1.5, n = 3 and |a]^ = A; the dashed curve is for the statistical 
mixture. 

In general for the superposition of coherent thermal optical fields in the cavity 
there are large scale macroscopic oscillations in the photon-number distribution which 
disappear for n >> \af, see |177j . In Fig. 5.11 we have plotted P(n) given by l |5.2.10| 
for TSEDS against the photon number n for fixed values of n, a and r (solid curve). 
Comparing our results with those obtained in |177j , we can say that P(n) in the present 
case have oscillations more pronounced than those given in the thermal even coherent 
state case. This emphasizes the role of squeezing parameter in the thermal optical 
cavity. In fact these oscillations are not only indication to nonclassical effects, but are 
also emerging from the quantum interference between the components of the state in 
the phase space. This can be seen when comparing dashed curve (corresponding to 
the statistical mixture) and solid curve in the figure (these curves are well defined for 
n positive integer or zero). 

Now let us discuss the degree of purity Trp^ for TSDS. Substituting l |5.2.4t into 
l |4.3.6| , carrying out the integration and putting |e| = 1, we get 



Tr(p|) = ^,{2 -h exp (-^) -h 4exp [-e^^aH2n + 1 + 2^)] 

(5.2.14) 

X cos((/)) + exp [-4e^''a^(2n + 1)] cos(2(/))}. 
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As we mentioned before for a pure state Tr(p^) = 1 and for a mixed state Tr(p^) < 1. 
The authors of [74] have shown that the degree of purity of the input chaotic light states 
is unchanged by subsequent squeezing or displacement process. Further, in |177| it 
has been shown that the superposition of coherent thermal fields in a cavity is quite 
complicated, thermal odd coherent state is more chaotic than the thermal light, but 
for thermal even coherent state it was found that < Tvip'^g^) < 1. In comparing 
equation (29) in |177| with expression l |5.2.14| , we can say that the degree of purity 
for thermal odd (even) coherent state is larger than that of TSODS (TSEDS) owing to 
squeezing parameter r. In Fig. 5.12 we have plotted Tr(p|) for TSVDS against a for 
fixed value of n = 1 and different values of r (r = 0,0.6). From this figure it is easy 
to see that the inequality < Tr(p|) < 1 is satisfied as for thermal even coherent 
state, where we find that Tr(p|) starts to increase its value comparing to the value for 
thermal light until it reaches the maximum, then it starts again to decrease until it 
takes the value of thermal light, but this is reached after long period of a. Also we 
noticed that as r increases, the Tr(p|) rapidly decreases to the value for chaotic light. 




0.32 1 

0.00 1.00 2.00 3.00 

Alpha 

Figure 5.12: Trp^ for TSVDS for n = 1, r = (solid curve) and r = 0.6 (dashed curve). 

5.2.5 Phase property 

As is well known, the nonclassical properties of electromagnetic waves are progress- 
ively destroyed by presence of noise and losses. Therefore the effect of thermal noise 
on quantum phase measurement, especially of nonclassical light, has attracted only a 
little attention |75]. Here we consider this effect based on the quasiprobability distri- 
bution function technique, which is more convenient for this problem. So inserting 
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W(]8, s) given by l |5.2.6|l into equation l |4.5.9|l and using the identity 



exp(-ax^ - bx - c)dx = -^a/ — [1 - erf( ^_ )], 
2\l a 2^/a 



where erf (.) is the Gauss error function defined by 



erf(x) = ^ / exp(-y^)dy, 
Jo 



after minor algebra we arrive at 



(5.2.15) 



(5.2.16) 



P(^)(0- 



Ay 



rVA+i^ + Ecos(20) 



2v^ ^^^^ 



+ 



[1 - |e|2 + (1 + |e|2)erf(h2yhr)] 



(5.2.17) 



+ |e|exp(T3) 



cos.j)exp()^) 1^3^^ 



+ ^^^^^[sin(/) + ierf (ihsvTn) cos (/)] 



where 



hi = 



Ti = 



A+M + Ecos(2e) 



(AH 



5)2 -£2 



T2 



h2 = 



a[A+M+E] COS0 



(AH 



5)2-E2 cos(20)' 



A+i^+Ecos(2e)' 



T3 



C2c 



_C2sin0_ 
(A+i^)+E' 



- G. 



■-2-' 



(5.2.18) 



The phase distribution functions obtained by means of quasidistribution functions 
may possess some difficulties, such as the singularity using the P-function and the 
negative values using W-function for some states. However, these difficulties can be 
avoided by considering the case s = -1, corresponding to the Q-function which is 
always positive. 

The phase distribution for superposition of quantum states, say even and/or odd 
coherent states, has two-peak structure [179] . For the present case we have displayed 
in Figs. 5.13a-c the phase distribution for TSYDS. In Fig. 5.13a we plotted the phase 
distribution for the Yurke-Stoler state with r = 0, and n = 0. In this figure we can see 
two-peak structure for even (odd) coherent state, which has been replaced by single 
peak at = and two wings as — > +t{. This behaviour has been realized in studying 
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(c) 

Figure 5.13: Phase distribution for TSVDS for a) P{9,a) for r = 0,ri = 0; b) P(0,r) for n = 0,a = 1; c) 
P{0,r) for n = 2,a = 1. 

the odd binomial states, see |180j . The squeezing effect in the Vurke-Stoler optical field 
in the cavity, i.e. squeezed Vurke-Stoler state (n = 0), has also been displayed in Fig. 
5.13b. It is clear in our case that the bifurcation phenomenon appears in range of 6 
shorter than that for the squeezed state, where the central peak is at 9 = and the 
lateral peaks are at = +1.5. Finally we have plotted in Fig. 5.13c the effect of thermal 
noise in addition to squeezing for fixed value of the parameter a (a = 2). From this 
figure we can realize that the latter bifurcation has been destroyed, e.g. the central 
peak is washed out, which indicates the influence of losses in the optical cavity, when 
the nonclassical effect disappeared. 

5.2.6 Conclusions 

In this section we have considered the thermal noise in the distribution of superposition 
squeezed and displaced number states. An observation of such a system can be found 
in a trapped ions, where a motion of a harmonically bound ^Be~^ ions can be used, 
see for example |35]. This can be seen when ions are initially cooled by a laser close 
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to zero point of motion and the coupling between their motional and internal states 
due to applied (classical) radiation can be described by the Jaynes-Cummings-type 
interaction. In this case the evolution of the internal atomic state will create states 
having quantum character, such as number state. On the other hand, if the ions are 
not in the Fock state, the motional state is characterized by a density operator whose 
diagonal elements have a number distribution P(n) given by 



where Pj is the probability of occupation in | |) (| |) is the hyperfine ground state of 
the ion), Qn,n+i is the Rabi frequency and Vn is the decoherence rate between levels 
|n) and \n + 1). P(n) in expression j5.2.19| has been derived analytically for more 
general case, see equation ( |5.2.10| , which can be controlled by the detuning Doppler 
cooling. Therefore we have an advantage to measure and control further phenomena 
by using the photon-number distribution given by ( |5.2.19| . 

For thermal superposition of squeezed and displaced states, the nonclassical effects 
have been traced via quasiprobability distribution functions, second-order correlation 
function, photon-number distribution, and phase distribution. We have seen that the 
correlation between different oscillators is essentially responsible for nonclassical ef- 
fects, which has been demonstrated by the behaviour of Wigner function and photon- 
number distribution. The thermal light photons have tendency to bunch each other, 
this has been reflected in the behaviour of our state, for example the normalized cor- 
relation function g'^^O) shows sub-Poissonian statistics but only for real and negative 
values of e, provided h is sufficiently small. Q-function exhibits always two peaks stru- 
cure provided that e f 0. Nevertheless, for the sake of comparison, it is important to 
remind that for the state jS.l.l^ . as we have shown in section 5.1, the W-function ex- 
hibits negative values for a wide range of phase space, i.e. the role of the interference 
term is less effective; g^^^O) shows sub-Poissonian statistics for a different type of e, i.e. 
e = +l,i; further Q-function displays multi-peak structures related to the excitation 
number n. 

We have also calculated the photon-number distribution using the Wigner function 
and we have shown that the oscillations in P(n) are increased when the squeezing is 
considered for superposition of thermal coherent fields in the cavity. Also we have 
found that the purity coefficient for this superposition which is mixed rather than pure. 
Finally, we have considered the phase properties of the optical system and shown that 
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for a short range of 6, bifurcation distinguishing the squeezed states can occur in the 
absence of thermal noise. This phenomenon disappears with increasing thermal noise 
in the cavity. 
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5.3 Quantum statistics of a solvable three-boson squeeze op- 
erator model 

In this section we introduce a highly correlated multidimensional squeeze operator. 
This operator consists of three modes in interaction and it represents the time-dependent 
evolution operator of the interaction part of the Hamiltonian |181II182| determined as 

f = Y!,]=i (^jajaj - iAi {did2exp[i(wi + co2)t] - h.c.} 

(5.3.1) 

-iA2 {did3exp[i((i>i + (i>3)f] - h.c.} - iX^ - d2d3 exp[i((i»2 - (1*3)^] - h.c. 
where dy and d| satisfy the commutation relations 

[ak,aj] = 6kj, Skj = (l l^lj, (5.3.2) 

and coj are the field frequencies and Xj are the effective intermodal coupling constants. 
In fact, Hamiltonian equation l |5.3.H describes a two-photon parametric coupling of 
modes 1 and 2, and 1 and 3 (two photons are simultaneously created or annihilated 
in both the quantum modes through the interaction with classical pumping mode), 
and linear interaction of modes 2 and 3. By introducing the transformation Ay = 
Aj exp{i(jOjt), the Hamiltonian l |5.3.1| may be transformed into the interaction picture. 



- = -iMiAiM - aIaI) - iA2(AiA3 - A[A3^) - iX^{AlA2 - A3A2^ (5.3.3) 

It is important mentioning that the interaction l |5.3.3t can be established in a bulk 
nonlinear crystal exhibiting the second-order nonlinear properties in which three dy- 
namical modes of frequencies (i>i, (i»2/ ti's are induced by three beams from lasers of 
these frequencies. When pumping this crystal by means of the corresponding strong 
coherent pump beams, as indicated in the Hamiltonian, we can approximately fulfil 
the phase-matching conditions for the corresponding processes, in particular if the 
frequencies are close each other (biaxial crystals may be helpful in such an arrange- 
ment). Also a possible use of quasi-phase matching may help in the realization, which 
is, however, more difficult technologically jl83j . Another possibility to realize such 
interaction is to use a nonlinear directional coupler which is composed of two optical 
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waveguides fabricated from nonlinear material described by the quadratic susceptib- 
ility x'^). Modes 1 and 2 propagate in the first waveguide and 1 and 3 in the second 
waveguide (auxiliary device such as bandgap quantum coupler [129j or a set of mirrors 
can be used to generate two identical modes (mode 1) in each waveguide). The inter- 
actions between the modes in the same waveguide are established by strong pump 
coherent light. In this case the coupling constants Ai and A2 are proportional to the 
second order susceptibility x'^^ of the medium and they also include the amplitudes 
of the pump. The linear coupling between modes 2 and 3 is established through the 
evanescent waves [96] . 

For completeness the time-evolution operator of equation l |5.3.3| is 

exp(-i^f) = exp[Aif(A[A; - A1A2) + A2f(A[A; - A1A3) + X^t{A^Al - AjAg)], (5.3.4) 
which can be identified well with a time-dependent three-mode squeeze operator 

S(r) = exp[ri (A[AJ - A1A2) + r2(A[A5 - A1A3) + r^iA^Al - AlAz)], (5.3.5) 

where = Xjt, with < rj < 00, j = 1,2,3 and r = (ri,r2,r3). It is evident that this 
squeeze operator must involve two different squeezing mechanisms (terms involving 
ri and rz) and then it is more complicated than squeezing operators that have appeared 
in the literature earlier ||2]|7]|39]- [431. Now if we set 

A = (A[AJ - A1A2), B = (A[A5 - A1A3), c = (A2A5 - A3AI), (5.3.6) 
then we have the following commutation relations 

[A,B] = -C, [B,C] = A, [C,A] = B. (5.3.7) 

We may conclude that the squeeze operator l |5.3.5t involves correlations and can be 
regarded as the exponential of linear combination of three generators, which are 
closed under the commutation relations j5.3.7|l , and it represents the su(l, 1) general- 
ized coherent state. As we mentioned earlier the quantum correlation between different 
quantum mechanical systems is responsible for the nonclassical effects for such type 
of operator. 
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Our plan of studying operator ( |5.5.5|l will be as follows: We derive the basic relations 
for three-mode squeeze operator and then we discuss three-mode photon-number 
sum and photon-number difference as well as squeezing phenomenon. Moreover, we 
give two examples of three modes squeezed states related to this operator which are 
three-mode squeezed coherent and number states. For these two states we discuss 
sub-Poissonian phenomenon and quasidistribution functions. 

5.3.1 Properties of the correlated quantum systems 

Squeezed state of light is distinguishable by long-axis variance of noise ellipse for one 
of its quadratues in phase-space. This property is connected with the pairwise nature of 
the unitary operator l |2.1.3| under which the initial state evolves. This operator exhibits 
some well-known basic relations summarized in the literature j2]. We introduce similar 
relations corresponding to the three-mode squeeze operator j5.3.5| , and then we use 
them to deduce the three-mode photon-number sum and difference. Also we discuss 
squeezing phenomenon related to this operator. 

The squeeze operator l |5.3.5|l provides a Bogoliubov transformation of the annihil- 
ation and creation operators that mixe the three modes as 

Ai s S~'{r)AiS{r) = Ai/i + A2V2 + Kh. 
As ^ S~\r)A2S{r) = Aagi + A^g2 + A\g^, (5.3.8) 
A3 s S-^{r)AiS{r_) = Ashi -h A[h2 + Mh^, 

where 

h = cosh;i + ||sinh2(H), = asinh^i - 2^sinh%). 
h = ^sinh;i + ^smh\^). g, = cosh^i 

g2 = ^sinh;i -h ^O^sinh^d), g^ = ^sinh;! + ^smh^^), (5.3.9) 
hi = cosh M - || sinh^d), hz = f sinh m - ^ sinh2(f ), 
h3 = ^sinh/-h^sinh2(f), 

where = + - r§ and r^ < rf + r^. The derivation of these relations is given 
in appendix B. 

The commutation relations l |5.3.2| under the transformations j5.3.8| hold also for 
the operators Aj. Using these transformations we can easily calculate the statistical 
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properties for each mode. It is worth to mention that the corresponding expressions 
for the single mode squeezed operator |2] and for two-mode squeezed operator ||39] 
can be obtained from l |5.5.8|l by taking r2 = = 0, and A3 — together with A2 — > Ai 
for only single mode case. 

We may point out that a strong correlation is built up between the three modes de- 
scribed by squeeze operator ( |5.3.5| . This is quite obvious for the case of the parametric 
amplification when two mode waves are mixed to generate a third wave via nonlin- 
ear medium, e.g. in an optical crystal with nonlinear second order susceptibility |2]. 
This can be demonstrated with the help of three-mode pure squeezed vacuum states 
■S(ll) 0/^=1 10;), where S(r) is the squeeze operator l |5.3.5| . In this case the eigenstates 
of the three-mode photon-number difference A[Ai - A2A2 - A3A3 correspond to zero 
eigenvalue, thus 

A(A[Ai - AJA2 - AjAs)^ = 0. (5.3.10) 

However, the situation will be different for three-mode photon-number sum, using 
l |5.3.8|l and after minor calculations we obtain 



A(A[Ai -h A^Aa + A^A3)2 = ff{ff - 1) + g|(i + g|) 

+ hi(l + hi) + 2(/fg| + ffhl + higf). ^ • ■ ^ 

In order to see the correlation between modes we have to calculate the fluctuations 
for both sum and difference operators thus having 



2 [{a%aIa2) + {aIa,a%) + {aIa^aIa,) - (A[Ai)(A5A2) ^ 
- (A[Ai)(A5A3) - (AJA2){AJA3)) = 2{ffgl + ffh^ + hjgl). 

It is clear that this quantity has non-zero value and this is the signature of photon- 
number correlation between modes. 

Now we show how quantum correlations between the systems can give rise to 
squeezing if operators act in the spaces not only corresponding to three systems (three- 
mode squeezing) but also to two systems (two-mode squeezing), rather than to the 
individual systems (single-mode squeezing). This will be done using three-mode pure 
squeezed vacuum states S(r) HLi |0/) as before. For this purpose, using results of 
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section 4.2, we calculate the squeezing variances for three-mode squeezing quadratures 
given by ( |4.2.4| and l |4.2.5| in terms of three-mode pure squeezed vacuum states. After 
minor algebra we get 



{(AX)2) = I [/2 + fl + fl + gl + gl + g^ 



(5.3.13) 



((Ay)2) = 1 [fl + fl + fl + gl + gl + gl 
+hf + hl + h^ + 4h2g3 - 4/ih2 - 4/ig3 



(5.3.14) 



As we mentioned earlier the expressions for the single-mode and two-mode squeez- 
ing can be obtained easily from l |5.3.13| and l |5.3.14| . It should be taken into account that 
C = 5/1/2 corresponding to the single-mode, two-mode and three-mode squeezing, 
respectively. 

From l |5.3.13t and l |5.3.14t one can easily prove that the model cannot exhibit single- 
mode squeezing, e.g. for the first mode we have 



Si = Ss = 2(/| -h fl), (5.3.15) 

where the relations between the coefficients fj resulting from the commutation rules 
of Aj have been used to get such relation. 

For compound modes the system can provide two-mode (only between modes (1, 
2) and (1,3)) as well as three-mode squeezing in the V-quadrature as shown in Fig. 5.14 
(for the shown values of squeeze parameters). From this figure one can observe that 
the behaviour of three-mode squeezing factors is smoothed in such a way that they 
are initially squeezed and their squeezing values reach their maximum, then they start 
again to decrease into unsqueezed values for large domain of r\ (which is not shown 
in the figure). Also, the initial values of squeezing are sensitive to the values of the 
squeeze parameters (compare solid and short-dashed curves). Concerning two-mode 
squeezing factor (long-dashed curve) one can see it is monotonically decreasing func- 
tion with lower limit -1. It is important mentioning that the squeeze operator under 
discussion cannot practically provide maximum squeezing (for three-mode squeez- 
ing), i.e. S2 = -1, and this, of course, is in contrast with the single-mode and two-mode 
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Figure 5.14: Sg against squeeze parameter ri for three-mode squeezing and for = 0.1, rg = 0.1 (solid 
curve) and 0.2 (short-dahed curves); and for two-mode squeezing (between modes 1 and 2, long-daslied 
curve) (r2,r3) = (0.2,0.1). 



squeeze operators ( |2.1.3| and ( |2.1.4| where they can display maximum squeezing for 
large values of squeeze parameter, however, they cannot provide squeezing initially. 
In conclusion it is quite obvious that squeezing can occur in the combined systems 
even if the individual systems are not themselves squeezed. The mechanism for this 
process is the correlation between the systems. 

5.3.2 Three-mode squeezed coherent and number states 

Here we examine the Glauber second-order correlation function for squeezed coherent 
states and squeezed number states related to the three-mode squeeze operator l |5.3.5| . 
Further we discuss the anticorrelation properties for three-mode squeezed coherent 
states. 

First we define the three-mode squeezed coherent states as 

|a,r) ^ S(r)D(a)|0)i|0)2|0)3, (5.3.16) 

where S(r) is the three-mode squeeze operator l |5.3.5t and D(a) is the three-mode 
Glauber displacement operator given by 

3 

D(a) = exp ^(ayA| - a* A,-), (5.3.17) 

7 = 1 
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and a = {a\,a2, a-^). 

The unitarity of the operator l |5.3.5[l provides that three-mode squeezed coherent 
states are not orthonormal but they are complete in the following sense 

i = j j j \a,r){r,a\d^aid^a2d'^a5. (5.3.18) 

Using l |5.3.8|l the mean photon numbers for various modes in three-mode squeezed 
coherent states are given by 



and the photon-number variances are 



((Ani)2),„h = (2/2 - l)(ni)coh - iff - 1)^ 
((An2)2)coh = (2g| + l)(n2)coh - §1 

where coh stands for squeezed coherent states. Expressions related with the 3rd mode 
can be obtained from those of the 2nd mode by using the following transformation 



{91.92.9^) (h3,hi,h2). (5.3.21) 

Having obtained equations ( |5.3.19| and ( |5.3.20| , we are ready to examine the second- 
order correlation function given by ( |4.1.H . In the following we restrict our discussions 
to the first mode 1, because the other modes would have similar behaviour. 

In phase space, squeezed coherent states \a,r) are represented by a noise ellipse 
with the origin at a, and they do not exhibit sub-Poisson distribution, i.e. they exhibit 
Poisson distribution at r = which is growing rapidly to superthermal distribution, i.e. 
g'2)(0) > 2, and it persists for a large domain of r f70ll71] . In our model of three-mode 
squeezed coherent states we show that they exhibit only partial coherence behaviour, 
i.e. 1 < g(2)(o) < 2. 

The condition for sub-Poissonian statistics is that the variance {{Ahj)^) must be less 
than the mean photon number (hj). We show that this condition will not be fulfilled 
for all modes, i.e. sub-Poissonian light cannot be obtained, because we have, e.g. for 
the first mode, the inequality 
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2\ffa, + hal + half + fl + fl < 0, 



(5.3.22) 



which will not be satisfied for any values of the coherent amplitudes aj. 
For the first mode Ai we have 



grio) = 1 + 



2{f! - l)(ni)coh - {f! - if 

("l)coh 



(5.3.23) 



It is clear that when Vj = 0, we recover the normalized second-order correlation 
function for coherent states. According to l |5.3.23| , to reach Poissonian statistics or 
thermal statistics, the expectation value (ni) of the photon number should have the 
following values: 



where subscripts p and th denote the corresponding quantities for Poisson and thermal 
distributions, respectively. For instance, to obtain thermal field, from ( |5.3.19| and 
( |5.3.24| , i.e. (ni)th' we have 



It is clear that l |5.3.25| is satisfied only when ay = 0, i.e. for three-mode squeezed 
vacuum states, and hence super-thermal statistics cannot be obtained. Similar pro- 
cedures show that Poissonian statistics cannot be obtained. This is a consequence 
of intermodal correlations of three-mode squeezed coherent states, for aj j= 0, light 
exhibits only partially coherence. This can be seen in Fig. 5.15, where we display 
the normalized normal second-order correlation function for the first mode against 
ri for aj = 2exp(i^),7 = 1,2,3, rs = 0.2 with r2 = 0.4,0.8, and 1.5 corresponding to 
solid, short-dashed, and long-dashed curves, respectively. Also we have displayed the 
corresponding normalized second-order correlation function for two-mode squeezed 
coherent state (centered curve) for the amplitudes aj = 2exp(i^),7 = 1,2, for sake of 
comparison. In this figure it is clear that partial coherence is dominant, and it persists 
for large values of ri and the initial values of gj (0) are sensitive to r2. However, the 
normalized second-order correlation function for two-mode squeezed coherent state 
(centered curve) is a growing function, starting from 1 (Poisson distribution) when 
there is no squeezing ri = 0, and becomes stable for large values of r displaying 



(ni)p = ^(/f - 1), (ni)th = /f-l. 



(5.3.24) 



l/itti + f2a*2 + half = 0. 



(5.3.25) 
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partial coherence behaviours. These values exceed those for three-mode squeezed 
coherent state for large domain of r. 




Figure 5.15: Normalized normal second-order correlation function gf^{0) (first mode) for three-mode 
squeezed coherent state against ri, when a,- = 2exp{i j),;' = 1,2,3, rj = 0.2, taking = 0.4 (solid curve), 
0.8 (short-dashed curve) and 1.5 (long-dashed curve); centered curve is corresponding to the second-order 
correlation function for two-mode squeezed coherent state and a, = 2exp(tj),j = 1,2. 



Second, we discuss the sub-Poisson properties of three-mode squeezed number 
states. This state can be written with the aid of three-mode squeeze operator as 



\n,r) s S{r)\ni)\n2)\n^), 



(5.3.26) 



for simplicity we set n = (ni, n2, 113). 

The mean photon numbers in three-mode squeezed number state are 



ni)n = niff + ("2 + 1)/! + ("3 + 
2 , ^ „2 



"2 



n2gf + nsgl + (ni + 



(5.3.27) 



and the variances of the photon number are 
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((Ani)2)„ = /f/|[ni + ns + 2nin2 + 1] + /f/|[ni + ns + Sni/lg + 1] 
+f2fi["-2 + ns + 2n2n3], 

(5.3.28) 

((An2)^)n = gfgiii^T, + "2 + 2n2n3] + gfgiifii +112 + 2nin2] 
+g|g|[ni + ns + 2nin3 + 1], 

where hj is the mean photon number for the 7th mode. The expression for the third 
mode can be obtained using l |5.3.21|l . It is worthwhile to refer to |70ll71|[73lll84lll85| , 
where further discussions related to squeezed number states are given. For instance, 
when the squeezing is not significant, i.e. r is small, the normalized second-order 
correlation function can be less than unity, which indicates that the light field has sub- 
Poissonian statistics |71]. Furthermore, squeezed vacuum state exhibits super-Poisson 
statistics (precisely superthermal statistics) for r f (24). Here, in contrast to the 
latter results, we prove that three-mode squeezed number state can exhibit thermal 
statistics when hj = 0,; = 1, 2, 3, i.e. for three-mode squeezed vacuum states, otherwise 
sub-Poissonian statistics or partially coherence behaviour is dominant. We focus our 
attention to the first mode. According to the definition of normalized second-order 
correlation function to obtain thermal statistics, we need to have 

((Ani)2)„ - (ni)„ - (ni)2 = 0; (5.3.29) 
from j5.3.27| and l |5.3.28| and after minor algebra one finds 

ni(ni + l)/i^ -I- n2(n2 + 1)/1 + n^{n^ + 1)/| = 0. (5.3.30) 

It is clear that the equality sign is only satisfied for fij = 0,; = 1,2,3. In other words, 
superthermal statistics will never occur. This is in contrast with the squeezed number 
state. 

In Fig. 5.16a we depict normalized second-order correlation function for single 
mode case for three-mode squeezed number state against r2 for fij = 1,; = 1,2,3 and 
(ri,r3) = (0.5,0.3). We notice in general that partially coherence behaviour is dominant 
again. Further we can see sub-Poissonian behaviour for small values of for modes 
1 and 3 with maximum value at the third one. For sake of comparison, we displayed 
the normalized second-order correlation functions for single mode (dashed curve) and 
two-mode (solid curve) squeezed number states in Fig. 5.16b for fij = 1,;' = 1,2, where 
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Figure 5.16: Normalized normal second-order correlation function for: a) three-mode squeezed number 
states against for different modes, when hj = 1,;' = 1,2,3, {ri.r^) = (0.5,0.3) for mode 1 (solid curve), 
mode 2 (short-dashed curve) and mode 3 (long-dashed curve); b) for single mode squeezed number 
states (solid curve, ni = 1) and two-mode squeezed number state (dashed curve, n,- = 1,/ = 1,2) against 
squeezing parameter r. 



we can see that both of them have sub-Poissonian statistics for lower values of squeeze 
parameter, otherwise they characterize partially coherence behaviour of light beams. 

Finally we turn our attention to discuss the effect of intermodal correlations in terms 
of anticorrelations in three-mode squeezed coherent states. We do this by investigating 
the behaviour of the Cauchy-Schwarz inequality giving by j4.1.6| . In doing so we have 
to derive the expectation values of cross photon-number operators between various 
modes, which are given in appendix C. In Fig. 5.17 we have plotted the quantity Jyj? 
indicating the violation of Cauchy-Schwarz inequality between the ;th mode and the 
kth mode against squeeze parameter ri, where (r2,r3) = (0.4,0.2) and a, = 1,; = 1,2,3 
(real). Further, solid, short-dashed, and long-dashed curves correspond to the Jyjj 
quantity obtained between (1, 2), (2, 3), and (1,3) modes, respectively. In general, photon 
antibunching can occur in dependence on the values of the parameters and a;. 
More precisely, in this figure we can observe that all curves can take on negative 
values reflecting the violation of the inequality, i.e. the photons are more strongly 
correlated than it is possible classically. The strongest violation of this inequality occurs 
in the (1,3) mode for lower ri and then the curve monotonically increases to positive 
values for larger ri. The weakest violation is in the (2,3) mode for which the curve 
decreases from positive values to the negative values as ri increases. Clearly, the 
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violation of this inequality is sensitive to the values of squeeze parameters and coherent 
amplitudes. As is known, the correspondence between quantum and classical theories 
can be established via Glauber-Sudarshan P-representation. But P-representation does 
not have all the properties of a classical distribution function, especially for quantum 
fields. So the violation of the Cauchy-Schwarz inequality provides explicit evidence of 
the quantum nature of intermodal correlation between modes which imply that the 
P-distribution function possesses strong quantum properties |150j . 

-t-> 0.10- 

D 

CT 0.05-, 




ri 

Figure 5.17: The quantity indicating the violation of the Cauchy-Schwarz inequality betwen the 
/th mode and feth mode, against squeeze parameter ri, where {r^.r^) = (0.4,0.2) and a,- = 1,/ = 1,2,3 
(real). Solid, short-dashed, and long-dashed curve correspond to the above quantity considered for modes 
(1, 2), (2, 3), and (1,3), respectively. 



5.3.3 Quasiprobabilify distribution functions 

Here we extend our work to include the quasiprobabilify distribution functions. We 
consider such functions for three-mode squeezed coherent and number states. We 
give only one example in this thesis for the joint quasiprobabilify functions for three- 
mode squeezed coherent states since this type of functions carries information on the 
intermodal correlations between different modes in the quantum mechanical system. 
Results of section 4.3 should be used. 

(a) Three-mode squeezed coherent states 
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Here we study quasiprobability distribution functions for three-mode squeezed co- 
herent state, specified by the density matrix Pcoh = \oc,r){lL>^\ ^i^d \^>L) is given by 
j5.3.16|l . Thus substituting Pcoh into l |4.5.1|l and using transformations l |5.3.8|l together 
with the Baker-Hausdorff identity for normal ordering of operators and after some 
minor calculations we have the s-parametrized joint characteristic function 



cfj,(Ci,C2,C3,s) = exp[i E;=i(s|C;f - IVjf)] 



(5.3.31) 



X exp[(a^?7i - air^l) + {a^r]2 - 0277^) + (a^r/s - 0377^)], 
where we have used the following abbreviations 

r?i = Ci/i - C2f73 - Clh2, r?2 = C2gi - CJ2 + C^h^. = C^h^ + C292 - Cih- (5.3.32) 

Having obtained the characteristic function, we are ready to find the s-parametrized 
quasiprobability functions for three-mode squeezed coherent state by inserting l |5.3.31|l 
into (|4.3.2t and evaluating the integral by using the identity (A.4) in appendix A three 
times; the calculations are straightforward but rather lengthy, thus we get 



Wfj,(j8i,i82,i83,s) 



TT^ DlD2(2/| + 2/| + l-s) 



■ exp 



'2/1+2/1+1-S 



X exp 



1 

4di 



(A2 - C2) cos^ (/) -h (B2 - D^) sin^ </) + (CD - AB) sin(2(/))l j (5.3.33) 



X exp 



1 

4V2 



(B2 - D2) cos2 (/) + (A2 - C2) sin^ 4) - {CD - AB) sin(2</)) 



where 



vi = p.1 cos^ 4> + P^2 sin^ </> + AI3 sin(2(/)), V2 = cos^ + Mi sin^ 4>-P-'5 sm{24>); (5.3.34) 
here 

;ii = i(2/| + 2/1 + 1 - s)-i [2/2(1 - s) - 2g|(l -h s) - (1 - s^)] , 

Ai2 = |(2/| + 2/f + 1 - s)-i [2/2(1 - s) - 2h|(l + s) - (1 - s2)] , (5.3.35) 

M3 = (2/1 + 2/2 -h 1 - s)-i [2/2g2hi + gih3(l + s)] . 

In equation ( |5.3.33| we have also used the following definitions 



Doctoral Thesis 



75 



Faisal Aly Aly El-Orany 



Chapter 5. Scientific results and their analysis: part I (static regime) 



A = (a* - as) + 2/1 g3 (2/1 + 2/| + 1 - s)~^{aX - ai), 
B = (a* - as) + 2hh2{2fl + If I + 1 - s)-i(a* - ai), 
C = (a* + as) - 2/ig3(2/| + 2/| + 1 - s)~\a\ + ai), 
D = (a* + as) - 2/ih2(2/| + 2/| + 1 - s)-i(a* + ai). 



(5.3.36) 



where <^ = 5 tan ( 7r=u~ )' '^z ^ ~ J^/')' ^ 1'2, 3 and aj represents the expectation 



value of the operators Aj given by equation (|5.3.8|l in terms of the coherent state. The 
three-mode functions W^qj^(j8i,]82,j83, s) given in l |5.3.33t are 6-dimensional Gaussian 
functions and display the nonclassical correlation nature by involving the terms A^, B^, 
etc. when the cross-terms J3i]82,j8ij83, etc. are not zero. Furthermore, we can see that 
the three-mode P representation does not exist at least for some values of squeezing 
parameters rj, i.e. if = or xig = or both (see equation j5.3.33| ), the physical 
reason for non-conditional breaking down of the P-function lies in the extremely 
strong correlation between the amplitudes of the three modes of the system during 
the evolution determined by squeeze operator f5]. Such tightly correlation causes 
that the modes may no longer fluctuate independently even in small amount such as 
allowed in a pure coherent states f59] . 

Now we turn our attention to single mode case, using similar procedure as before. 
The s-parametrized quasiprobability function takes the form 



Light fields for which the P-representation is not a well-behaved distribution will exhibit 
nonclassical features. Clearly single mode P-function is well defined, because the para- 
metric systems evolution broadens P-distribution (increases the radius of the Wigner 
contour compared to the initial one) |39] and reflects that there is no (single mode) 
nonclassical behaviour, e.g. sub-Poissonian statistics and squeezing. Furthermore, it 
has been shown that P-function for the superposition of two fields is the convolution of 
the P-function for each field considered individually |139] , so that j5.3.37| (with s = 1) 
describes the superposition of P-function of a coherent state with complex amplitude 




(5.3.37) 



where 7 = 1,2,3 and 



Ti = 2/1 -1- 2/1 -hi, T2 = 2g| + 1, T3 = 2h| + 1. 



(5.3.38) 
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ocj and P-function of a chaotic mixture with variance ^{tj - 1) jS], i.e. displaced thermal 
light. Further, equation (|5.5.37t has a Gaussian form in phase space with width and 
center dependent on rj, aj and with a circular symmetric contour as a result of the 
fact that two quadrature variances are equal. Consequently the single mode photon- 
number distribution does not exhibit oscillations, in contrast with this for squeezed 
coherent state |168j , owing to the noise ellipse of W-function which is isotropic. 

Now we investigate single mode phase distribution in terms of Q-function, i.e. s = 
-1. Using the same procedures as in section 5.2.5 we arrive at 



P(0) 



+ 



1 + erf(- 



exp 

V2bj 



2(t; + 1) 



V(t; + l)exp 



2(t^ + 1) 



(5.3.39) 



where 



bj = Uj exp{-i9) + ajexp(i0), (5.3.40) 

with 7 = 1,2,3, ocj have the same meaning as before, and erf(.) is the Gauss error 
function l |5.2.16|l . 




Figure 5.18: Phase distribution P(0, ri) against 9 and ri for the first mode as output from three-mode 
squeezed coherent state for a,- = 1,; = 1,2, 3, ^2 = 1 and = 0.5. 

Here we restrict our investigation to real ay, i.e. by = 2ay cos 6, detailed examination 
to the formula l |5.3.39| shows that it is a 2;T-periodic function, a symmetric function 
(P(-0) = P(0)) around = and also it has its maximum height at 9 = 0. Moreover, 
this formula has a similar structure as that for coherent states (which can be recovered 
from l |5.3.39| by setting ry = 0). This means that the phase distribution of the single 
mode exhibits a one-peak structure for all values of ay and ry (see Fig. 5.18 for shown 
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values of parameters). That is the phase distribution of the single-mode as output from 
three-mode squeezed coherent state is insensitive to the quantum correlations between 
the systems. This situation is the same as that of two-mode squeezed coherent state, 
i.e. when rg = rs = |61]. However, it has been shown that the joint phase distribution 
for the two-mode squeezed vacuum depends only on the sum of the phases of the two 
modes, and that the sum of the two phases is locked to a certain value as the squeezing 
parameter increases |61ll62j . 

(b) Three-mode squeezed number state 

Here we calculate single mode s-parametrized characteristic and quasiprobability func- 
tions for mode 1, for three-mode squeezed number state l |5.3.26|l . Therefore s-parametrized 
single mode characteristic function can be written as follows 

ci^{Cu,s) = exp[i|Ci|2(s -2/1 -2/|)]L„,(|Ci|Vf)L„,(|Ci|2/|)L„3(|Ci|2/|), (5.3.41) 

where l |5.1.4|l has been used, Lk{.) are the Laguerre polynomials of order k, and sub- 
script sn stands for three-mode squeezed number states. 

The single mode s-parametrized quasiprobability function for mode 1 is given by 
Fourier transformation of l |5.3.4H as 

where ti is given in l |5.3.38| . In derivation of l |5.3.42|l we used differentiation under 
the sign of integral to perform the integration and the identity (A.4) in appendix A. 
Expressions for the second mode and third mode can be obtained from l |5.3.42| by 
replacing functions by functions gj and hj. 

Here we discuss the behaviour of W-function for the first mode when it is in the 
Fock state |1) while the other modes are in vacuum. In fact, W-function of squeezed 
Fock state 1 1 ) is well characterized by inverted hole which is stretched as a consequence 
of squeezing parameter [7lJ. In the following we can use Li(x) = 1 - x to analyze the 
behaviour of the model under discussion. For this case it is clear that, from equation 
( |5.3.42| , W-function can exhibit negative values only inside the circle < -^p- 

with center at the origin. However, the maximum value will be established at the 




(5.3.42) 
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Figure 5.19: W-function for the single mode (mode 1) as output from three-mode squeezed number 
state, assuming the first mode is |1) and the other modes are in vacuum for {ri.ri.r^) = (0.8,0.9,0.6). 

circle |j8ip = -^j^- Hence, W-function will not exhibit stretching since the variables 
RejSi and Im]8i are not involving squeezing factor which plays an essential role for 
stretched quadratures of squeezed number state. As a result the radii of these circles 
are dependent on squeeze parameters rj, so the negative values of W-function will be 
sensitive to the values of squeeze parameters, i.e. the function has negative values for 
a range of |]8i| shorter than for individual Fock state where it exhibits negative 
values inside the circle < ^. This can be recognized if one compares the well 
known shape of W-function of Fock state |1) with that of the single-mode case as output 
from three-mode squeezed number states. For this purpose Fig. 5.19 is displayed for 
W-function of |1,0,0) against x = Re]8i and y = ImjSi for (ri,r2,r3) = (0.8,0.9,0.6). 
We can see that the well known negative values of W-function of Fock state |1) are 
smoothed out, in turn, the circle of negative values is enlarged, as we have shown 
earlier. 

The phase distribution for the ;th mode, which is in Fock state \nj) while the other 
modes are in vacuum states, evolves under the action of three-mode squeezed operator 
j5.3.5| , with the aid of Q-function, to 




(5.3.43) 



which is uniform distribution. In other words, in spite of the system is highly cor- 
related the phase distribution in the single mode of purely nonclassical state is phase 
insensitive. This agrees with the fact that the states are represented by a density matrix 
which is diagonal in the number state basis, having random phase distribution |58] . 
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5.3.4 Conclusions 

In this section we have introduced new type of multidimensional squeeze operator 
which is more general than usually used and which includes two different squeezing 
mechanisms. This operator arises from the time-dependent evolution operator for 
the Hamiltonian representing mutual interaction between three different modes of the 
field. The origin of the nonclassical effects of this operator model is the correlation 
between the systems where we have shown that the quadratures squeezing can occur 
in the combination systems rather than in the individual systems. 

The quantum statistical properties corresponding to this operator have been traced 
by means of the variances of the photon-number sum and difference, squeezing phe- 
nomenon, Glauber second-order correlation function, violation of Cauchy-Schwarz 
inequality, quasiprobability distribution functions and phase probability distribution, 
considered for three-mode coherent and number states. 

For three-mode squeezed coherent state, we found that the second-order correla- 
tion function describes partially coherent field provided that one mechanism of squeez- 
ing is always surviving, which can be demonstrated also by means of quasiprobability 
distribution function in single modes. Nevertheless, this behaviour is in contrast with 
behaviour of squeezed coherent states, where second-order correlation function can 
display superthermal statistics, i.e. g'^^O) > 2. We have found strong violation for 
Cauchy-Schwarz inequality in some modes, i.e. the photons are more strongly cor- 
related than it is allowed classically. Concerning the single-mode phase distribution, 
single peak structure is dominant for all values of parameters provided that coher- 
ent amplitudes are real. The signature of the correlations between the three modes 
appears straightforwardly in the form of quasiprobability functions. 

For three-mode squeezed number states, the second-order correlation function is 
in agreement with that for the squeezed number state and in general it exhibits partial 
coherence; however, sub-Poissonian behaviour is attained for small values of p. and 
the maximum value is obtained only for = 0,; = 1,2,3. Nevertheless, it cannot 
exhibits superthermal statistics. Also the range of negative values of the W-function, 
for a single mode, are highly sensitive to squeeze parameters. Phase distribution, for 
a single mode, is insensitive to correlation between modes, and it displays a uniform 
form. 
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In this part we investigate the quantum properties of some dynamical systems in the 
framework of interaction of radiation modes with nonlinear media described nonlin- 
ear couplers. We have considered three types of such models when two, three or 
four modes are evolved within. In general, all of these models are composed from 
two waveguides of the length L through which the interacting modes evolve. The 
interaction between the different modes in the same waveguide can be established 
by strong laser pumping, however, the interaction between the two waveguides oc- 
curs by means of evanescent waves. Instead of initial coherent states at the input 
beams, the Fock states generated for instance in a micromaser can be introduced. 
Outgoing fields are detected as single or compound modes by means of homodyne, 
photocounting or coincidence detection. Further, the values of the coupling constants 
between the different modes in the device are chosen for the reason of simple calcula- 
tions and illustration of the physical behaviour of the system. We mainly demonstrate 
the effect of their mutual relations and their magnitudes are not substantial for such 
demonstrations. For an experimental realization these values can be estimated as 
(10^^ - 10^^)s~^, provided that the pumping laser of power in mW is used, producing 
about 10^^ photons/s, which is sufficient power to neglect quantum noise in pumping 
beams. Then all the effects shown are interpreted on the corresponding reduced time 
scale. Initial state amplitudes or photon numbers are chosen in units because switching 
properties are examined on the quantum level of single photons. 



81 



Chapter 6. Scientific results and their analysis: part II (dynamic regime) 



In the following section we investigate the quantum properties of the first model 
which includes two-mode case. 



6.1 Quantum statistics and djznamics of nonlinear couplers 
with nonlinear exchange 

In this section we discuss the quantum statistical properties of nonlinear couplers 
composed of two waveguides operating by the second subharmonic generation as- 
suming strong coherent pumping and linear exchange of energy between waveguides 
by means of evanescent waves, however we additionally take into account the influence 
of nonlinear coupling of the parametric type of both the waveguides. This system is 
suggested to be described by the following Hamiltonian: 



J- = X]/=i ■ (^i&jcLi + aj'^ exp{i^jt) + h.c 



+A3 



did2exp[i(/)i(0] + h.c. | + A4 {aid2 exp[-i(j)2{t)] + h.c.} , 



(6.1.1) 



where di(d[), a.2{a.l) are annihilation (creation) operators of the fundamental modes 
in the first and second waveguides having frequencies coi and co^, respectively, p.j are 
related with the frequency of the second-harmonic modes described classically as 
strong coherent fields, (f)j{t), j = 1,2, are related to the difference- and sum-frequencies 
of modes 1 and 2, respectively, Ai and A2 are nonlinear coupling constants for the 
second subharmonic generation in the first and second waveguides, respectively, A3 
is the coupling constant for linear exchange between waveguides through evanescent 
waves, A4 is the coupling constant for the nonlinear exchange through simultaneous 
annihilation or creation of a photon in both the subharmonic modes on expense of 
pumping and h.c. means the Hermitian conjugate terms (for further details concerning 
the optical parametric processes, see ^8) (Chap. 10)). 

When = and only the degenerate term is considered, we have the well-known 
Hamiltonian, in the interaction picture, for squeezed light generation [7], where Ai 
(or A2) represents the coupling constant proportional to the quadratic susceptibility 
of the second-order nonlinear process (degenerate parametric down-conversion with 
classical coherent pumping), or the coupling constant proportional to the cubic sus- 
ceptibility of the third-order nonlinear process (degenerate four-wave mixing with 
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classical coherent pumping) |186] . If additionally (^i(0 = (j)2W = 0- l^he Hamiltonian 
( |6.1.H represents a mixture of second subharmonic generation, frequency conversion 
and parametric amplification in the interaction picture P ISTlllOO] , Schematically, this 
Hamiltonian is represented in Fig. 6.1. 
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Figure 6.1: Scheme of quantum nonlinear coupler with linear and nonlinear coupling 
formed from two nonlinear waveguides described by the quadratic susceptibility x'^^- 
The beams are described by the photon annihilation operators as indicated; z = vt is 
the interaction length. Both the waveguides are pumped by strong classical coherent 
waves. Outgoing fields are examined as single or compound modes by means of 
homodyne detection to observe squeezing of vacuum fluctuations, or by means a set of 
photodetectors to measure photon antibunchibng and sub-Poissonian photon statistics 
in the standard ways. 

Solution of the equations of motion in the Heisenberg picture for the Hamiltonian 
(|6.1.1t are 

ai{t)exp{-itl^) = di(0)Ki(f) + df(0)Li(f) + d2(0)Mi(f) + dJ(0)JVi(f), 

(6.1.2) 

d2(0exp(-iff ) = d2(0)K2(0 + dJ(0)L2(f) + di(0)M2(f) + df(0)N2(f), 

where the explicit forms of the time dependent coefficients, which contain all the 
features of the structure, and the details about the solution are given in appendix D. 

One can see from this solution (see appendix D) that when and Q2 are real, 
the coupler switches the energy between the modes which propagate inside since the 
solution will include trigonometric functions [5]. Nevertheless, if Qi and Q2 are pure 
imaginary, the Heisenberg solutions provide hyperbolic functions, which are growing 
rapidly with time, and the coupler operates as amplifier for the input modes |187] . So 
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that the behaviour of the coupler will be indicated essentially by the relation between 
coupling constants. 

In what follows, we shall employ the results obtained here to treat the squeezing 
phenomenon and quasiprobability distribution functions for the model under consid- 
eration. 

6.1.1 Squeezing phenomenon 

Here we study single-mode squeezing phenomenon when the modes are initially pre- 
pared in thermal states with the average thermal photon numbers hj, j = 1, 2 as well 
as in the coherent states. For this purpose (using results of section 4.2 and considering 
lij{t) to be the phase of the local oscillator in the relations l |4.2.4|l and j4.2.5| , without 
loss of generality, we can cancel the high frequency terms, and ; = 1,2 stands for 
mode 1 and mode 2, respectively) we calculate the quantities Sj{t) and Qj{t) which are 
defined in (4.9) for the initial thermal light. For the first mode we have the following 
expressions 

Siit) = 2ni[|Li(f)|2 + |Ki(f)|2] + 2n2[mt)f + m{t)f] + + 2|Ni(f)|2 

(6.1.3) 

+2(2ni -h l)[Li(f)Ki(f) + c.c] + 2(2n2 + l)[Mi(f)JVi(f) -h c.c], 

Qi(f) = 2ni[|Li(f)|2 -h |Ki(f)|2] + 2n2[|Ni(f)|2 + \M,{t)f] + 2|Li(f)|2 -h 2\Ni{t)f 

(6.1.4) 

-2(2ni + l)[Li(f)Ki(f) + c.c] - 2(2n2 + l)[Mi(f)JVi(f) + c.c], 

where c.c. means the complex conjugate terms. The corresponding expressions for 
the second mode can be obtained from l |6.1.3| and l |6.1.4| by using the interchange 
1 2. Further, the other expressions related to the injected coherent light in the 
coupler are the same as j6.1.3| and l |6.1.4| , but setting hj = 0. 

It is known that the nonlinear coupler is a source of optical fields, the statistical 
properties of which are changed as a result of the linear and nonlinear interaction 
inside and between waveguides. Consequently, one can obtain nonclassical light from 
one input and, in addition, it can be switched. 

We have plotted Si(f), Qi{t) in Figs. 6.2a,b and Szit), Q2{t) in Figs. 6.3a,b, when the 
initial light is coherent, for different values of X^. Further we have chosen A3 = 1 for all 
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a) 



b) 




Time 



Time 



(b) 



Figure 6.2: Squeezing phenomenon for mode 1 when the modes are initially in coherent light and: a) 
for the first component Si(f); b) for the second component Qi{t); A3 = 1 for all curves; curves A, B and 
C are corresponding to Ai = A2 = A/, = 0.25, Ai = Ag = A4 = 0.2 and Ai = 0.17, A2 = A4 = 0.2, respectively; 
straight line has been put to show the bound of squeezing. 



curves and for the curve A : Ai = A2 = A4 = 0.25; for the curve B : Ai = A2 = A4 = 0.20, 
and for the curve C : Ai = 0.17, A2 = A4 = 0.2. On the other hand. Fig. 6.4 gives Si(f) 
(first mode) when the initial light is thermal light with coupling constants as those for 
the curve C, where h\ = 0.5 and fi2 = 0.5 (solid curve), 1.5 (dashed curve); straight line 
shows the bound of squeezing of the curves. Firstly, we start our discussion by studying 
the case of input coherent light. From these figures we can see how the coherent states, 
which are minimum-uncertainty states, evolve in the coupler to produce squeezed light. 
We can observe the oscillatory behaviour in these curves, showing that squeezing can 
be switched from one waveguide to the other in the course of time during power 
transfer. Moreover, squeezing can be interchanged between the two quadratures of 
the same waveguide. More precisely, for mode 1, squeezing can occur for all selected 
values of A^ in S\{t), but in Qi(f) only curves A, B can exhibit squeezing, as shown in 
Figs. 6.2a,b, which reflects the dependence of nonclassical behaviour on the strength 
of subharmonic generation. 

For mode 2, we can see squeezing in all curves in both the quadratures, as shown 
in Figs. 6.3a,b. It can be easily seen that the amount of squeezing is sensitive to the 
strength of coupling Ajj and that in general its values in the second component are 
more pronounced than those in the first one. Now if we turn our attention to the 
case of injected thermal light (Fig. 6.4c), we can observe that squeezing is available in 
the large interaction time. Further, Si(f) exhibits oscillatory behaviour and it evolves 
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from unsqueezed values in the short range of interaction time, owing to the fact that 
thermal states are not minimum-uncertainty states, into squeezed values and eventually 
unsqueezed values can be recovered. Indeed, we proved numerically that this beha- 
viour is periodically recovered with the time. Moreover, by comparing the dashed 
curve with the solid one, we can see that increasing of the photon number in the 
second waveguide causes decreasing of the amount of squeezing in the first one. This 
is related with the effect of evanescent waves between waveguides and shows how one 
can control light by light in the coupler. Finally, we can conclude that by controlling 
the input average thermal photon number and the interaction time (or the length of 
the coupler), the interaction under consideration can generate squeezed thermal light. 
It is worthwhile to refere to |70II71II75|[74] , where more discussions related to squeezed 
thermal states are given. 



6.1.2 Quasiprobabilify functions 

Here we shall continue in our investigation for the statistical properties of the system 
under discussion in the basis of quasiprobability distribution functions for single-mode 
case only when both the modes are initially in number and coherent states. 

(a) Input Fock states 
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2.40 - ;' 
2.00 - ;' 
1.60 - ; 




Time 

Figure 6.4: Squeezing phenomenon for mode 1 when the modes are initially in thermal light for the 
first component Si(t) with hi = 0.5, hi = 0.5 (solid curve), 1.5 (dashed curve) and the coupling constants 
Xj are the same as those for the curve C in Fig. 6.2. 

The density operator for two-mode number states is 

Pn(0) = |n)i|m)22(m|i(n|, (6.1.5) 

where the subscript n stands for Fock states case. The single-mode s-parametrized 
characteristic function for the first mode can be obtained using the results of section 
4.3 together with the identity ( |5.1.4| as 

d^UCsJ) = exp |||C|' - ^[\Vi{t)f + \V2{t)f]'^Ln{\vM%m{\V2{t)f), (6.1.6) 

where 

viit) = cm) - cLiW. V2{t) = cMr(f) - cm). (6.1.7) 

and Lm(.) is the Laguerre polynomial of order m. Inserting l |6.1.6t into l |4.5.4|l . carrying 
out the integration numerically (the exact form can be found, see jl89] ) and taking 
s = we obtain the W-function. The result has been plotted in 

Figs. 6.5 against x = Re]8i and y = ImjSi, when the first state is the Fock state |1) and 
the second one is the vacuum state |0), i.e. n = l,m = 0; A3 = l,Ai = A2 = A4 = 0.25 
and for shown values of time. In Fig. 6.5a we have the W-function for t = i.e. 
after short time interaction between the two modes we observe similar behaviour as 
for the W-function of the state |1), which means that pronounced negative values are 
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a) ^ b) 




(b) 

c) 




(c) 



Figure 6.5: W-function for the single mode (mode 1) for different values of time f when both the modes 
are initially in the number states; the first mode is in the state ]1) and the second mode is in the state ]0) 
and Aj = 1, Ai = Ag = A4 = 0.25: a) for f = 7r/100; b) for t = 7x12; c) for t = jr. 

exhibited. This behaviour of the W-function is completely different by increasing the 
time (f = ^); we see disappearance of negative values of the quasidistribution and a 
stretched positive peak occurs (Fig. 6.5b). This form of W-function is close to that of 
squeezed vacuum states [7], i.e. squeezed vacuum states can be generated, in principle, 
in this model. It should be borne in mind that the specific direction of stretching 
for the quasiprobability function of squeezed states may be achieved by choosing a 
suitable value for the phase of squeeze parameter. Of course, in Fig. 6.5b, there is 
a negligible spike at the top of the peak which can be smoothed out by governing 
the coupler parameters. After larger time interaction t = tt, the negative values are 
reached again, but they are less pronounced and asymmetry can be observed due to 
stretching (Fig. 6.5c). So we meet a time development of the W-function as a result 
of the power transfer between the two modes inside the coupler. This shows that the 
nonlinear directional coupler can be used as a source of quantum states. Further, we 
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have found that the behaviour of the Q-function is consistent with that of W-function 
for the same shown values |189] . 

(b) Input coherent light 

In a similar way as we followed in the case (a) we can study the same quantities when 
both the modes are initially in coherent states |ai, ag). For the single-mode case the s- 
parametrized characteristic function and the s-parametrized quasiprobability function 
are given, respectively, as 

cfi(Cs, t) = exp {-i[l - s + 2|Li(f)|2 + 2|Ni(f)|2]|C|2 + Cat(f) - C*ai(f)} 

X exp{^emmm + mm)]} (6.1.8) 

X exp {iC*^[JVi(f)Mi(f) + Li(f)Ki(f)]} , 



X exp 



X exp 



7rV[¥ + ii'lWP + lNi(f)[2]2-|Bi(f)l2 

[l^ + |Li(0|g + [Ni(f)p]|ai(f)-^r 
[i^ + lLi(f)l2 + |Ni(f)|2]2-[Bi(f)|2 



(6.1.9) 



^ + |Ll(f)[2 + [Nl(f)[2]2-|Bi(f)]2 



where 



Bi(0 = Nl{t)M*{t) + Limiit) = |Bi(f)|exp[2i<5i(f)], 
Ei(f) = [di(f) -J8]exp[i6i(f)], 



(6.1.10) 



and ai(f) is the average of the first mode operator Ai(f) with respect to input coherent 
light. 

From l |6.1.9|l . the P-function is not well defined as an ordinary function for \Li{t)f + 
|JVi(f)|^ < |Bi(f)| and hence the nonclassical effects such as squeezing of vacuum fluc- 
tuations can occur, as we have seen before. 

Furthermore, the nonclassical effect, especially squeezing of vacuum fluctuations 
in the case of our system, can be recognized in the behaviour of W-function (and/or 
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a) b) 




(a) (b) 

Figure 6.6: W-function for the single mode (mode 1) for different values of time t when both the modes 
are initially in the coherent states; |ai|^ = laal^ = 2 and Aj, are the same as in Fig. 6.6: a) for f = jt; b) for 

t = 27T. 

Q-function) in phase space as shown in Fig. 6.6 for shown values of parameters. For 
f = 0, i.e. when there is no interaction between the two modes, the W-function is 
identical with that shown for a single mode representing a symmetric Gaussian bell 
in phase space. As soon as the interaction switches on [t > 0), we observe that the 
Gaussian centre is shifted and the rotationally symmetric function of the initial state 
at f = gets to be squeezed in various phase space directions in dependence on time, 
as demonstrated in Figs. 6.6a,b. In other words, the initial symmetric contour of the 
W-function has been stretched as the interaction switches on, i.e. noise ellipse char- 
acterizing squeezed light appears, which rotates in the phase space as the interaction 
time progresses. 

6.1.3 Conclusions 

In this section we have examined the quantum statistical properties of radiation gen- 
erated and propagated in the nonlinear optical coupler composed of two nonlinear 
waveguides operating by the second subharmonic processes, coupled linearly by evan- 
escent waves and nonlinearly by nondegenerate optical parametric process. We have 
demonstrated regimes for generation and propagation of nonclassical light exhibited 
by squeeezing of vacuum fluctuations. We have also obtained quasidistribution func- 
tions for the initial light beams which are in Fock and coherent states. Compared to 
earlier results for nonlinear optical couplers we have shown that the nonlinear coup- 
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ling increases in general quantum noise in the device even if in some cases it can 
support generation of nonclassical light. 

The motivation for examination of the system under consideration arises from the 
previous investigations of the nonlinear couplers as promising devices to produce non- 
classical light. When coherent light is injected initially in the system, squeezed as well 
as sub-Poissonian light can be generated |ll89j (the latter has not been analyzed here). 
For injected number states, squeezed vacuum states are produced. When thermal light 
initially enters the coupler, the coupler can operate as a microwave Josephson-junction 
parametric amplifier |63]. These effects have been recognized to result from the com- 
petition between linear and nonlinear properties of the system and are dependent on 
the initial amplitudes of the input fields. The crucial role plays here the mechanism 
of the energy exchange between waveguids. 
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6.2 Generation of squeezed light in a nonlinear asszmmetric 
directional coupler 

In this section we study the properties of three modes interacting in a nonlinear asym- 
metric directional coupler composed of a linear waveguide and a nonlinear waveguide 
operating by nondegenerate parametric amplification. We show that such a device 
is a source of single-mode squeezed light. This fact has been demonstrated, under 
certain conditions and for specific modes, for incident coherent beams in terms of the 
quasiprobability functions, photon-number distribution and phase distribution. 

6.2.1 Equations of motion 

We consider a device (nonlinear asymmetric directional coupler) as outlined in Fig. 
6.7. It can be described by the Hamiltonian 



^ = (Ajjajdj] + iXi {did2 exp[i((L»i + (i)2)t] - h.c.} 



+ iXo 



did3exp[i(a)i - o)z)t] - h.c. 



+ iA?5 



d2d3exp[i((i»2 - (i>3)f] - h.c. 



(6.2.1) 



where aj{aj), j = 1,2,3 are the annihilation (creation) operators designated to the 
signal, idler and linear modes, respectively, Ai and Xj, ; = 2, 3 are the corresponding 
nonlinear and linear coupling constants; coj are the natural frequencies of oscillations 
of the uncoupled modes and h.c. is the Hermitian conjugate. The linear exchange 
between the two waveguides establishes through evanescent waves provided that a>i - 

The dynamics of the system is described by the Heisenberg equations of motion 
which, using the slowly varying forms (dj = Ayexp(-ia»yf), ; = 1,2,3), read 



^ = -XiAl - A2A3, 

^ = -X,aI-X^A^, (6.2.2) 
^ = A2A1 + A3A2. 
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These are three equations with their Hermitian conjugates forming a closed system 
which can be easily solved, restricting ourselves to the case A2 = A3 = to avoid the 
complexity in the calculations; this means that we consider stronger nonlinearity in 
the first waveguide compared with the linear exchange between waveguides, thereby 
having the solution 



Ai(f) = di(0)/i(f) + df(0)/2(f) - d2(0)/3(f) 

-dj(0)/4(f) - d3(0)/5(0 - al{0)U{t), 



Mt) = d2(0)gi(f) + ai{Q)g2{t) - di (0)^3(0 
-d[(0)g4(f) - d3(0)g5(0 - &liO)g6{t), 



(6.2.3) 



A3(f) = d3(0)hi(f) + d3^(0)h2(f) + d2(0)h3(f) 

+dl{0)hi{t) + di(0)h5(f) + al{0)he{t), 

where the time-dependent coefficients, i.e. fj{t),gj{t),hj{t), including all information 
about the system, are 



d2(0) 



di(0) 



d3(0) 



X 



(2) 



Ai 



A2 



X 



(1) 



A3 



jCLii 



Figure 6.7: Scheme of realization of interaction in l |6.2.H using a nonlinear asymmetric 
directional coupler which is composed of two optical waveguides fabricated from first- 
order (x^^^) and second-order (x'^^) materials, where x designate susceptibility. Signal 
mode 1 and idler mode 2 propagate in the first waveguide and linear mode 3 in the 
second waveguide. The interaction between the signal and idler modes is established 
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by strong pump coherent light, which is not indicated in figure, with the coupling 
constant Ai. The interactions between the linear mode and the signal and idler are 
established linearly with the coupling constants A2 and A3, respectively. 



/i,3(0 = I cosh(Aif) + cosh(^)cos(fef) + ^ sinh(^) sin(fef) 



/2,4(^) = I Isinh(Aif) +sinh(^f)cos(fef) + ^ cosh(^) sin(fef) 



1 



/5(0 = vlcosh(^)sin(fef), /6(f) = 



/§sinh(^)sin(fef). 



(6.2.4) 



hi(f) = cosh(^) cos(fef) - ^sinh(^)sin(fef)' 
h2{t) = -sinh(^)cos(fef) + cosh(^) sin(fef). 



^t3,5(f) = COSh(^)sin(fef), h4,6(0 



/§sinh(^)sin(fef). 



where k = ^Ai and the expressions for gj{t) are the same as those of fj{t). 

In fact, the nature of the solution can show how the coupler does work. To be more 
specific, the time-dependent coefficients contain both trigonometric and hyperbolic 
functions. Consequently, the propagating beams inside the coupler can be amplified 
as well as switched between waveguides, i.e. between modes, in the course of time. 

On the basis of the well known commutation rules for boson operators, the follow- 
ing relations can be proved for the time-dependent coefficients 



ff it) - flit) + flit) - flit) + flit) -flit) = 1, fiii)g,it) - f2ii)gm 



+/3(0g2(f) -fdi)gi{i) - m)gdi) + HW) = o. 



fi[i)g^{i) - f2[i)gS) + h{i)9i{i) - h{i)92{i) - f^{i)g5{i) + f^{i)g^{i) = 0. 



(6.2.5) 



The remaining relations can be obtained from j6.2.5| by means of the following trans- 
formations 
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^ {-g^W. -9aW. -9iW. -g2W.g5{t).9Q{t)) (6.2.6) 
^ (^h5{t),he{t), -h^{t), -hi{t), -hi(f), -h2{t)). 

Based on these results, we can study the quantum properties of the evolution of 
different modes in the model when they are initially prepared in coherent states. This 
will be demonstrated in the following parts by means of quasiprobability functions. 

6.2.2 Quasiprobability functions 

In the following we consider phase space distributions for the single-mode when all 
modes are initially prepared in coherent states \ai,a2, as) before entering the coupler. 
We start from knowledge of the single-mode s-parametrized characteristic function 
which is suitable to describe the quantum statistics of light and for the first mode it 
takes the form 

Ci(Cs,f) = exp |i|C|2[s - |T/i(f)|] + |T?2(f)|(C' + C') + [«t(OC - ai(OC*]J . (6-2.7) 
where rij{t), j = 1,2, are given by 

= /f (0 + flit) + flit) + flit) + flit) + flit), 

(6.2.8) 

r?2(0 = flit)f2it) + fmdt) + f5it)fGit). 

Now the s-parametrized single-mode quasiprobability functions of the signal mode, 
which can be derived in a similar way as we did in the previous section, are 



Wi(j8,s,f) = 



X exp ■ 



jr^[A+{t)~s][A4t)~s] 



(6.2.9) 



2[A_(f)~s] 2[A+(f)-s] 

where li{t) = ai(f) - j8 (ai(f) has the same meaning as before) and 
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A±(f) = [hit) ± /2(f)f + [/3(f) ± hWf + m) ± /6(0f ; (6.2.10) 

the other expressions for the idler and linear modes can be obtained from ( |6.2.7| and 
( |6.2.9| by making use of the transformations ( |6.2.6| : we except here the case s = 1, 
which will be discussed shortly. Now let us start our investigation by demonstrating 
the behaviour of the W-function, i.e. setting s = in l |6.2.9| . In the language of 
mechanical analogy of a harmonic oscillator with dynamical conjugate variables x and 
y, i.e j8 = X + iy, the W-function can be written in the form 



W.(x,y,0^ , ■ exp ] - ^1:1^1 K (6.2.11) 



2^\/((AX(())2>((Af(()F> i 2((AX(())2> 2{(AVm) 
where 

((AX(f))2) = ((Ay(f))2) = (6.2.12) 

andai(f) = ay-{t)+iay{t); ((AX(f))^) and ((Ay(f))^) are single-mode x- and y-quadrature 
variances, respectively. 

Now from l |6.2.11t and l |6.2.12t together with l |6.2.4t one can easily prove that the 
quadrature variances of the linear mode (single-mode squeezing) are 

((AX(f))2) = i ||sin2(fef) -h [cos(fef) -h -^sin(fef)]2| exp(-Aif), (6.2.13) 

{{AV{t)f) = ^ 1 1 sin2(fef) + [cos{kt) - sm{kt)f | exp(Aif), (6.2.14) 

where k has the same meaning as before. Thus squeezing can be achieved in the 
linear waveguide in the X-quadrature. Moreover, it is clear that squeezing values 
become more pronounced for a large interaction time, i.e. for long length of the 
coupler L. This behaviour shows that changing the power of the linear interaction it 
is possible to transfer nonclassical properties, which are generated in the nonlinear 
waveguide, to the linear signal mode. For completeness, the uncertainty relation reads 

1 16 

((AX(f))2)((A V(f))2) = ^[1 + Y smHkt)]. (6.2.15) 
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This formula reveals that the minimum-uncertainty relation holds only when t = tg = 
where m is positive integer and we call tg a squeeze time. In this case the 
device provides the squeezed coherent light in the linear mode with squeeze para- 
meter r = It is important to mention that at i = ig the two waveguides become 
completely independent (cf. equations HSZ)) and the signal and idler modes can 
display perfect two-mode squeezing, which one can easily check with the help of two- 
mode quadrature variances. This can be explained as follows: A portion of the energy 
always remains within the guide into which the field was initially injected. This energy 
grows in the course of time until t = tg, when the two guides are completely independ- 
ent and hence the modes are trapped in their own guides. So that the linear mode 
can yield the well known single-mode squeezed light in the linear waveguide having 
its origin in the nonlinear waveguide in which the nonlinearity produces two-mode 
squeezing. For later times t j= tg, the device generates single-mode squeezed light 
in the linear waveguide where the corresponding W-function may be broader than 
that of the idealized squeezed coherent states [7]. This situation is close to that of a 
two-photon absorber (a two-photon absorption by a reservoir of two-level atoms from 
a single mode of the electromagnetic field jl90] ) where a squeezed state, which is not 
a minimum-uncertainty state, has been generated fl91] . The origin of squeezing of 
initially unsqueezed light interacting with two-photon absorbers is that the squeezing 
is generated by simple quantum superposition of states of light |190| . 

We proceed in our discussion by focusing our attention on the P-function. We 
have shown earlier that this system is able to provide squeezed light and this should 
be reflected in the behaviour of the P-representation, i.e. setting s = 1 in j6.2.9| . The 
significant example for this situation is the behaviour of the linear mode. For this case, 
we can show that the single-mode quadrature squeezing is established provided that 
A+(f) - 1 < and A_(f) - 1 > 0. It is evident that the P-function is not well-behaved 
function in this case and this is the indication of the nonclassical fields. Indeed, if we 
look at the model under consideration of a competition of parametric processes, we 
can find that this result is in contrast to the results in the literature for simpler systems 
of interacting modes starting from coherent states |(5l [59ll40ll45lll81] . where the delta 
function of the initial P-representation of coherent light for a single mode becomes 
well-behaved distribution during the interaction. 

In the following we use the phase space distribution functions to study the photon- 
number distribution and phase distribution for the system under discussion. 
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6.2.3 Photon-number distribution 

Photon-number distribution, i.e. the probability of finding n\ photons in the signal 
mode at time t, can be obtained in the photodetection process, and can be determined 
by means of the relation l |4.3.8| as 



P(ni,f) 



V[A+(f) + l][A-(f) + l] 



exp 



2[A_(f) + l] 



[ai{t]+al{t)f 
2[A+(f) + l] 



r=0 



A+(f) + l 



A-(0-l 



A„(f) + 1 



ni ~r 



(6.2.16) 



■(l+A+(f))(A+(f)-l) 



T "5 
^ni ~r 



(l+A„(f))(A-(f)-l) 



where LJJ^ (.) is the associated Laguerre polynomial and the integration has been carried 
out using the same technique as in appendix A. 

It is convenient to compare this distribution with the corresponding Poisson distri- 
bution 



P(ni,f) = ^ ' exp[-(ni(f))], (6.2.17) 
ni! 

which corresponds to fully coherent field with the same mean photon number {hi{t)). 
The expressions for the idler and linear modes should be obtained from j6.2.16| using 
the transformations ( |6.2.6| . 

In Fig. 6.8, we have plotted P(ni, t = 1.5) against ni for the signal mode (solid curve) 
for aj = 3exp(i2.) and Ai = 0.3. For the sake of comparison, the corresponding photon 
distribution for coherent field, (16.2.17^ . is shown by dashed curve. For all numerical 
calculations 5Zru=o ^ ('^i' = 1 with good accuracy. We noted for the signal mode that 
the photon-number distribution exhibits always one-peak structure. From this figure 
we see that the behaviour is rather super-Poissonian as a result of quantum fluctuations 
because the solid curve is always broader than the corresponding dashed curve. So 
that, in general the Poissonian light evolves in the nonlinear waveguide as super- 
Poissonian light. On the other hand, we have realized earlier from the behaviour of 
VV-function that squeezed states can be generated. However, these states are exhibiting 
oscillating photon-number distribution for certain values of squeeze parameter. 

Indeed, these oscillations can be recognized for the linear mode where this mode 
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25.00 



ni 

Figure 6.8: Photon-number distribution P(ni) for the signal mode for a,- = 3exp(iy),; = 1,2,3, t = 1.5 
and Ai = 0.3. The dashed curve is the Poisson photon-number distribution. 

can display squeezed light (see Fig. 6.9, for shown values of parameters). From 
this figure we see the macroscopic oscillations of the photon-number distribution for 
squeezed light. Further, comparison of Figs. 6.9a and 6.9b shows that the nonclassical 
oscillations of P(n3) are faster and more pronounced when t increases. It is worthwhile 
mentioning that the case of Fig. 6.9a corresponds to squeezed light with squeeze 
parameter r ~ 3.6. 

In fact, the behaviour of the photon-number distribution in this device is slightly 
different from that of the nonlinear asymmetric coupler with strong classical stimu- 
lating light in the second harmonic mode |120| where oscillatory behaviour as well 
as sub-Poissonian statistics for specific modes in both linear and nonlinear waveguide 
are exhibited. 

6.2.4 Phase distribution 

We make use of the single-mode Q-function to investigate the phase distribution for 
the system under discussion using (4.20). The result of the integration has the form 



P{9,t) = 



■C(f)VA+(f)A4f) 



1 



exp 



2[A_{()-H] 2[A+(f)-H] 



(6.2.18) 




exp[||l[l+erf(J|=j)l}, 
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b) 



35 40 AB 




(a) 



(b) 



Figure6.9: Photon-number distribution P(n3) for the linear mode when a, = 0.5exp(i^),/ = l,2,3;Ai = 1 
and a) t = ; h) t = 5. 



where 



C(f) 



r[A+(0 sin^ + A_(f) cos2 0], 



A_(t)A+(t)L--^^ ■ (6 219) 

^(^) = A„(f)l(f) [-^-(^) + - ^^+(^) - cos 0] , 

and erf (.) is the Gauss error function; A± (f) = 4((AX± {t)f) with ((AX(f))2) and ((A 

are single-mode x- and y-quadrature variances, respectively; di(f) is the expectation 

value of the signal mode operator in coherent state. 



i 



(b) 
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Cl, 



m 




(a) 
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Figure 6.10: Phase distribution P{0, t) for the signal mode against 6 and f for Ai = 0.6 and a) ay = 0; b) 
a, = 1, ;• = 1,2,3. 
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The behaviour of this distribution in this device can be understood on the basis of 
the competition between the two-peak structure for the vacuum states and the single- 
peak structure for coherent states provided that the coherent amplitudes are real. 
Such behaviour is demonstrated in Figs. 6a,b for the signal mode and in Figs. 7a,b 
for the linear mode, where the evolution of the phase distribution j6.2.18| has been 
depicted against the phase 9 and the time f. 




On the other hand, the well-known behaviour for the phase distribution of squeezed 
states may be established when we focus our attention on the behaviour of the linear 
mode, as expected (see Figs. 6.11a and b). From Fig. 6.11a one can see that the 
two-peak structure is typical for squeezed vacuum states, i.e. we have two-peak struc- 
ture for 9 = +^ IETIES]. Further, the heights of these peaks are ^\J^^ any 
time f (f > 0). However, Fig. 6.11b displays the well-known bifurcation shape for 
the phase distribution of squeezed field, i.e. the distribution curve undergoes a trans- 
ition from single- to a double-peaked form with increasing time. Indeed this figure 
qualitatively is quite similar to that for idealized squeezed coherent states fSO], and the 
distributions differ in the behaviour of the initial peak, which is here amplified for a 
while before splitting into two peaks. This is connected with the significant influence 
of the development by nonlinear effects and power transfers between waveguides. It 
is worthwhile mentioning that a similar behaviour has been obtained for a contra- 
directional nonlinear asymmetric coupler with strong stimulated coherent field in the 
second harmonic waveguide |150] . 
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Finally, we would like to conclude this part by discussing the origin of such beha- 
viour of the phase distribution. This can be easily understood by analyzing the function 
P{0, t) when aj = 0, ; = 1,2,3. In this case the formula l |6.2.18|l reduces to 



VA+{t)A4t) 
2;TA+(f)sin2 + A_(Ocos2 0' 

This distribution function is of double- or three-peak structure according to the relation 
between A+{t) and A_(f). Let us restrict our discussion to the one of the modes which 
are propagating in the nonlinear waveguide, say, to the signal mode. For this mode, it 
can easily be shown that A+{t) > A_(f) and consequently the formula l |6.2.20t exhibits 
three-peak structure with peaks for 6 = 0,+7r. The heights of these peaks at any time 
f f are ^ XITJ- non-zero displacements ap an additional factor of the form 
for coherent state, but with a coherent amplitude ai(f) which is the expectation value 
of the signal mode operator in coherent state, appears in the distribution. This factor 
is responsible for a peak at = related to coherent component, which competes 
with the three-peak structure of vacuum to provide the previous behaviour of the 
distribution. Similar argument can be adopted to explain the behaviour of the linear 
mode, however, it should be borne in mind that in this case A+(f) < A_(f), as we have 
shown earlier for squeezing phenomenon, and hence the two-peak structure for the 
input vacuum case is dominant at = 



6.2.5 Conclusion 



In this section we have shown that squeezed light can be generated in a nonlinear asym- 
metric directional coupler composed of a linear waveguide and a nonlinear waveguide 
operating by nondegenerate parametric amplification. After using the Heisenberg ap- 
proach to the quantum statistics of interacting modes, we have investigated such an 
effect in the linear mode, propagating in the linear waveguide, in terms of the qua- 
siprobability functions, photon-number distribution and phase distribution when the 
modes are initially in coherent states. Concerning the behaviour of Glauber P- and 
Wigner W-functions for the linear mode we can conclude that the former is not well- 
behaved function when f > 0, while, the latter displays the well-known behaviour for 
squeezed light in which one of the quadratures is amplified and the other is attenuated. 
For the photon-number distribution, the large scale macroscopic oscillations related to 
squeezed light are established. The phase distribution displays the bifurcation typical 
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for squeezed field. In fact, this behaviour shows that changing the power of the linear 
interaction, it is possible to transfer nonclassical properties, which are generated in 
the nonlinear waveguide, to the linear signal mode. Further, we have shown, in gen- 
eral, that the system has less tendency to generate sub-Poissonian light from initial 
coherent light. 
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6.3 Quantum sfafisfical properties of nondegenerate optical 
parametric S57mmetric coupler 

In this section we shall concentrate on studying the statistical properties of an op- 
tical field propagating within a directional coupler containing a parametric amplifying 
medium. Our starting point is the Hamiltonian, which represents a nonlinear direc- 
tional coupler composed of two nonlinear waveguides (taking into account the case of 
a strong pump, where non-depleting intensity is included into the amplifier coupling 
constant). This model can be described by an effective Hamiltonian as 

— = Ai(did3 + alal) + A2(did2 + d[d2) + A3(d2d4 + dgdj), (6.3.1) 

and the corresponding scheme is illustrated in Fig. 6.12. Two waveguides of the 
length L operating by nondegenerate optical parametric processes with signal beams 
described by annihilation operators di and d2 and idler beams described by ds and dj^ 
are pumped by strong laser beams of complex amplitudes aip and a2p, respectively and 
signals are coupled linearly by evanescent waves with the strength A2; v is the speed 
of the waves (dispersion is neglected) and is quadratic susceptibility. The coupling 
constants Ai and A3 include the pumping amplitudes. It is important to mention that 
in Fig. 6.12 the beam splitter plays only auxiliary role to incline a strong classical 
beam from a laser and consequently the vacuum fluctuations in the free port can be 
neglected. Both the beams before entering the waveguides must be attenuated to a 
low-intensity quantum level (the attenuators are not indicated in the figure) and only 
these beams are involved in the described quantum interaction. 

Clearly, the Hamiltonian l |6.3.1[l has a kind of symmetry which will be helpful to 
reduce the number of equations. For example, if one takes di < — > d2, ds < — > dj, and 
Ai < — > A3, the Hamiltonian l |6.3.1t is still invariant. 

The above Hamiltonian has been considered earlier by Janszky et al. j95] (taking 
into consideration the coupling parameters Ai and A2 are equal), where they discussed 
the propagation of a quantum field in coupler when the channels contain a parametric 
amplifying medium. This situation occurs if one or both channels |193 j are made from 
a second-order nonlinear material realizing a parametric process (down conversion). 
They have found two operation regimes: one under threshold, where the amplific- 
ation is less effective than the coupling between channels, and the other regime is 
above threshold, where the amplification constant is greater than the coupling con- 
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stant. However, the above Hamiltonian l |6.5.1|l will give us an advantage to see spatial 
effects (e.g. switching between the channels), where the coupling parameters will play 
a great role in controlling this phenomena. This will be demonstrated in the following. 
Also we shall give more details for the squeezing phenomenon and sub-Poissonian 
effect of this model [931. 



Signal laser 




A2 



z=vt 



Idler laser 



^a4( 



Figure 6.12: Scheme of realization of interaction in (6.31); BS are beam splitters. 



6.3.1 Equations of motion and their solutions 

The exact solution of the Heisenberg equation of motion for the bosonic operators of 
the model are 
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AiW = kimdO) + 1/2(^)02(0) - ih{t)al{0)+h{t)al{0). 



&2{t) = g2{t)&2{0) + igi(f)di(0) + ^3(^)03^(0) - igiimliO). 



(6.3.2) 



d3(f) = h^m^iO) + h2{t)d2{0) - ihi{t)a[{0) - ihi{t)a[{0), 



ai{t) = U(f)d4(0) + li{t)a[{0) - il2{t)a^{0) - ihWa^iO), 



where dy(0) are the input operators and the expression of the time dependent functions 
are given in the appendix E. 

From equations l |6.3.2| together with the exact forms for the time-dependent coeffi- 
cients in the appendix E, it is easy to check that the well known commutation relations 
between bosen operators are satisfied. 

In the following we treat some of the quantum statistical properties of the model 
under discussion, using the results presented here, related to squeezing phenomenon, 
sub-Poissonian and antibunching phenomena. 

6.3.2 Two-mode squeezing phenomenon 

Here we study two-mode squeezing using quadrature operators X and V which are 
defined in a similar sense as in section 4.2. Since the electromagnetic field is guided in- 
side the structure, exchange of energy between the two waveguides is possible because 
of the evanescent field between the waveguides [194] . So we can examine several cases 
of two-mode squeezing for this system using the input initial coherent states 0^=1 1";)- 

In the first case, two-mode squeezing for mode 1 and mode 4, i.e. for the signal 
mode in the first waveguide and the idler mode in the second waveguide, is described 



by 




sin(fQi) 
Hi 



2 



(6.3.3) 



2((AV)2)exp(-2(^) = [cosh(^cos(ff^i) - sinh(^cosh(ff^2)]^ 
+ f (A3 cosh (/) - A2 sinh 0) ^n^^M _ (A3 sinh (p - X2 cosh (p) ^ 



sin(fQi) 



2 



(6.3.4) 
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In the second case, two-mode squeezing for mode 1 and mode 2, i.e. for the signal 
modes in the first and second waveguides, is described by 

((AX)2) = {{AVf) = i[l + /32(f) + flit) + glit) + g|(f)]. (6.3.5) 

In the third case, two-mode squeezing for mode 3 and mode 4, i.e. for the idler 
modes in the two waveguides, is described by 

((AX)2) = {{AVf) = 1[1 + hi(f) -h hi(f) + [?(0 -h llit)]. (6.3.6) 




We can see from 1 16.3.3^ -1 (6.3.6^ that this coupler can generate squeezing, in terms of 
two-mode definition, only when considering the signal mode in the first waveguide and 
the idler mode in the second waveguide or vice versa. So we concentrate on this case. 
In Fig. 6.13, we have plotted squeezing component for the first case for different values 
of Xj against time t. Curves Al and A2 are related to the X-component corresponding 
to the two groups of values for Xj respectively, and similarly curves Bl and B2 are 
related to the V-component. For Ai small, A2 > A3 > Ai, i.e. when linear exchange 
between waveguides is stronger than nondegenerate parametric amplification inside 
waveguides, we observe that squeezing is dominant in the first quadrature having 
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oscillatory behaviour; nevertheless the value of squeezing is negligible in the second 
component, see curves Al and Bl. Increasing the values of both Ai and A2 and keeping 
A3 as before, we observed squeezing in both the quadratures and its values were more 
pronounced than in the earlier case with the maximum value of squeezing in the first 
quadrature, see curves A2 and B2. 

Thus we can conclude that this system is able to generate squeezed light in terms of 
two-mode squeezing, only when the signal and idler modes are considered in different 
waveguides, provided that the signal modes are coupled. The values of squeezing are 
well controllable by the values of coupling constants. 



6.3.3 Second-order correlation function 

As we have mentioned earlier sub-Poissonian photon statistics need not be associated 
with antibunching, but can be accompanied by bunching. We trace the nonclassical 
effects for the model under discussion using the normalized second-order correlation 
function for the various modes when the initial input light modes are in number 
states ftj^i well as in the coherent states. Then, we extend our discussion to 

demonstrate the photon antibunching phenomenon for our model using the relation 



(4.4). 



When the input light is in the number states, we have 



(ni(f + r))n = ff{t + T)ni + + r)n2 + fl{t + T)(n3 + 1) 
-h T)(n4 + 1), 



(6.3.7) 



(n3(f + T))n = hl{t + T)n3 + hf (f + T)n4 + h|(f + r){n2 + 1) 
+ hf(f + T)(ni + 1), 



(6.3.8) 



and 



Doctoral Thesis 



109 



Faisal Aly Aly El-Orany 



Chapter 6. Scientific results and their analysis: part II (dynamic regime) 



{almlit + r)a,{t + T)di(0)n = IffWn, + fl{t)n2][fl{t + r){n^ + 1) 
+/|(f + T)(n4 + 1)] + hMfmm + t) + h{t)f2{t + t)]2 
+/2(f)/2(f + T)ni(ni - 1) + fl{i)fl{i + T)n2(n2 + 1) 

+ 2[/i(f)/i(f + T)ni + /2(f)/2(f + T)n2][/3(f)/3(f + T)(n3 + 1) 
+/4(f)A(f + T)(n4 + 1)] 

+ [/2(f + T)ni + /|(f + T)n2][/|(f)(n3 + 1) + /|(f)(n4 + 1)] 

+ (n3 + l)(n4 + i)[m)fS + t) + /4(0/3(^ + t)]2 

+/|(f)/|(f + T)(n3 + l)(n3 + 2) + /|(f)/|(f + T)(n4 + i){n, + 2), 



(d3^(f)d3^(f + T)d3(f + T)d3(f))n = [hl{t)h^ + hl{t)nj,][hl{t + T)(n2 + 1) 

+hf(f + T)(ni + 1)] + n3n4[h3(f)h4(f + t) + hj,{t)h'i{t + t)]^ 
+h2(f)h§(f + T)n3(n3 - 1) + hf(f)h|(f + T)n4(n4 + 1) 

+2[h-i{t)h-i{t + T)n3 + hj,{t)hj,{t + T)n4] ^ ^ 

x[h2(f)h2(f + T)(ni + 1) + hi(f)hi(f + T)(ni + 1)] 

+ [hi(f + T)n3 + + T)n4][hf(0(ni + 1) + h2(f)(n2 + 1)] 

+ (n2 + l)(n2 + l)[/^l(^)/^2(^ + t) + h2(f)hi(f + r)f 

+hl{i)h\{i + T)(ni + l)(ni + 2) + hf (f)h|(f + T)(n2 + l)(n2 + 2), 

where the subscript n stands for the Fock states case. 

Firstly, we will discuss the sub-Poissonian statistics, i.e. we take t = (see the 
discussion in section 4.1) in the equations ( |6.3.7| -l |6.3.10| . For iij = 0, i.e. for input 
vacuum states, we can get 



((Ani(f))2)o - (ni(f))o = (/|(f) + /|(f))2, 
((An3(f))2)o - {Um = {h\{i) + hl{i)f, ^ • • ^ 



where the relation such as 



mhj{t)f) - {hj{i)) = (df (f)d2(f)) - {a]{i)aj{t)f, (6.3.12) 

has been used. It is clear from j6.3.1H that for input vacuum states, which have max- 
imum pronounced sub-Poissonian statistics, they evolve into pure super-Poissonian 
statistics states for all times f f inside the coupler. We display the second-order 
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Figure 6.14: Normalized normal second-order correlation function when both the modes are initially in 
the number states with mean photon numbers hi = A.hj = 1,/ = 2,3,-4, for all curves Aa = 1.2, A3 = 0.5 
and for curve A, Aj =0.1 and for curve B, Ai = 0.6: a) gf (f) for mode 1; b) grf*(t) for mode 2, straightline 
is corresponding to g'^^O) of ^he coherent light. 



correlation function (4.1) against time t in Fig. 6.14 for input number states with dif- 
ferent mean photon numbers hj {h\ = 4, n,- = 1,; = 2,3,4) and for different values of 
coupling constants. 

(2) 

In Fig. 6.14a, {t) for mode 1 has been shown, where always A2 = 1.2, A3 = 0.5 
and Ai = 0.1,0.6 are corresponding to curve A and curve B, respectively. Initially, 

(2) 

gl (0) = 0.75, which is the corresponding value for the Fock state |4). Increasing 

(2) 

the time, switching between modes starts (see curve A), g^ (f) has an oscillating be- 
haviour between values for sub-Poissonian and super-Poissonian statistics for rather 
short interaction times. However, for large interaction times the oscillations are suc- 
cessively washed out and super-Poissonian statistics are dominant. As seen in curve 

(2) 

B, gl '{t) does no exhibit any periodic behaviour and takes values corresponding to 
super-Poissonian statistics shortly after the switching on of the interaction. On the 
other hand, we can see that strongly sub-Poissonian light is generated not only for 
short time interaction but also for larger times as demonstrated by the behaviour of 

(2) 

g2 {t) in the second mode for the same values of rij as in the former case and for 

(9) 

(Ai,A2,A3) = (0.1,1.2,0.5), see Fig. 6.14b. Moreover, g^ \'^) behaves more smoothly 

f21 

than g[ (f) for the same values of hj and Ay. 

We note that for A2 less than Ai and A3, f f the second-order correlation func- 
tions exhibit always super-Poissonian statistics and there is no oscillatory behaviour 



Doctoral Thesis 



111 



Faisal Aly Aly El-Orany 



Chapter 6. Scientific results and their analysis: part II (dynamic regime) 



regardless of the values of input mean photon numbers, since all time dependent func- 
tions in equations (E.1)-(E.2) of appendix E are being hyperbolic functions, which are 
monotonically increasing functions. 

a) 



T3 





00 10.00 20,00 30.00 40.00 50.00 60.00 0.00 10.00 20.00 30.00 40.00 50.00 60.00 

Time Time 



(b) 



Figure 6.15: Normalized normal second-order correlation function gj (f) for different modes when both 
the modes are initially in the coherent states with |a,] = 1,/ = 1, 2, 3, 4, = f and ipj = ^.j = 2,3,4 for 
all curves, Aj = 0.1,A2 = 1.2, A3 = 0.5: a) first mode-curve A and fourth mode-D; b) second mode-curve B 
and third mode-C. 



For input coherent light, we have 

(ni(f + T))coh = \alfi{t + t) - ia^/2(f + t) + ia^h{t + t) + aj,fi{t + r) 
+fl{t + r)+f!{t + T), 



(6.3.13) 



{h^{t + T))coh = |a3^i3(^ + t) - ialhi{t + t) + alh2{t + r) - iaj,hj,{t + r)\^ (53 14) 
+hf{t + t) + h2(f + t). 



(d[(f)d[(f + T)di{t + T)di(f))coh = (ni(f))coh(ni(f + T))coh 
Hhmii + r) + hmii + r)f 

+ MW + r)+ h{t)h{t + T)][(df(f))coh(ai(f + T))coh + c.c]. 



(6.3.15) 



{al{t)al{t + T)a^{t + T)d3(f))coh = (n3(0)coh(n3(f + T))coh 

-h[hi(f)hi(f + T) + h2(f)h2(f-hT)]2 

+ [hi{t)hi{t + t) + h2(f)h2(f + T)][(d5(f))coh(a3(^ + T))coh + c.c]. 



(6.3.16) 
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where (ny(O)coh can be obtained from {hj{t + T))coh by simply setting t = 0; (dj(f))coh 
is the expectation value for d|(f) in the coherent states, c.c. is the complex conjugate; 
and coh stands for the expressions which are related to input coherent light. From 
equations (4.1) it is clear that the condition for sub-Poissonian statistics is that the 
variance {{Ahj{t))'^) is less that the mean photon number {hj{t)). Combination of 
( |6.3.13| - j6.3.16| into l |6.3.12| , after taking t = 0, it is easy to show that this system cannot 
provide sub-Poissonian light for input coherent light for t > 0. For instance, for mode 
1, from equations j6.3.12| with ( |6.3.13t and ( |6.3.15| it follows that sub-Poissonian light 
could be generated provided that 

2\alh{t) - ialh{t) + iaMi) + "4/4(0 1^ + /f (0 + /f (0 < 0. (6.3.17) 

It is evident that this inequality cannot be fulfilled regardless of the values of aj and 
Xj. More precisely, for input coherent light, this model can generate classical light, e.g. 
coherent, partially coherent and chaotic light. 



a) b) 




(a) (b) 



Figure 6.16: Normalized normal second-order correlation function gr] [t) for different modes when both 
the modes are initially in the coherent states with ]ai] = 10, = i.j = 2,3,4, rpj = = 1,2,3,4 where 
A; have the same values as in Fig. 6.15: a) first mode ; b) third mode-curve B and fourth mode-C. 

(2) 

In Figs. 6.15 and 6.16 we have plotted {t) against time t for Ai = 0.1, = 
1.2, A3 = 0.5 and input coherent light with complex amplitudes aj = \aj\e^i. In all these 
figures gf\o) = 1 initially. In Fig. 6.15 for |a;| = 1 and (ipi.ipz.^z-^i) = (f,f,f,f), 

(2) 

we see the oscillatory behaviour in the evolution of (f) in all cases showing that the 
photons are transferred from one mode to the other. These oscillations are success- 
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ively disappearing for large interaction time. They can be destroyed by changing the 
values of the complex amplitudes of input light, but not by changing the values of the 
coupling constants (see Fig. 6.16 for shown values). This shows how one can control 
light by light via nonlinear medium. In Fig. 6.16a we can see the periodical behaviour 
exhibiting photon statistics between close to Poissonian and super-Poissonian states 
for mode 1 with gradual increase of values of g^^HO = 2 representing chaotic light 
for larger interaction times. We see that the initial coherent state can be success- 
ively approximately regenerated. From Fig. 6.16b we can see that the precise Poisson 
distribution is obtained only initially and partially coherent light is obtained for later 
times, with the maximal noise value in mode 3 (curve B). For input coherent light 
2 > gy (f) > 1 holds, i.e. one cannot obtain superchaotic light. One can also observe 
some complementary behaviour of both the modes. 




-1.30 1 

0.00 5.00 10.00 15.00 20.00 25.00 

Time 

Figure 6.17: The quantity Ki(t) for mode 1 when both the modes are initially in the number (solid curve) 
and coherent (star-centered curve) input states for the same situation as in Fig. 6.14a (curve A) and Fig. 
6.15a for number state and coherent state, respectively. 

Concerning photon bunching and antibunching according to the definition (4.4), we 
have analysed the cases of input number and coherent states. We note, in general, that 
the quantity Kj{t) exhibits oscillatory behaviour between negative and positive values 
for both number and coherent input states, i.e. both antibunching and bunching can 
occur. The photon antibunching is more pronounced for input number states than for 
input coherent states. Comparing these results with those for sub-Poisson statistics, we 
can conclude that there is no direct relation between antibunching and sub-Poissonian 
statistics here, in agreement with results shown in the literature earlier |1 45111 44] . We 
can demonstrate this graphically in Fig. 6.17, by considering the first mode, for input 
number state (solid curve) and for input coherent state (star-centered curve) for the 
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same values of parameters as those of Figs. 6.14a (curve A) and 6.15a for number 
state and coherent state, respectively. 

6.3.4 Conclusions 

The main conclusions of this section can be summarized as follows. Two-mode squeez- 
ing, sub-Poissonian and antibunching effects of an optical field propagating inside a 
directional coupler containing nondegenerate parametric amplification process have 
been studied in the framework of Hamiltonian formalism. Incident number states and 
coherent states have been considered. We have shown that for input coherent light, i.e. 
Poissonian light, the system can generate squeezed light in terms of two-mode squeez- 
ing depending on the values of the coupling constants. More precisely, when the linear 
coupling is stronger than the nonlinear one, there is possibility to obtain squeezed light 
considering signal beam in one waveguide and the idler beam in the other, provided 
that the signal beams are connected by the evanescent waves. Nevertheless, this pos- 
sibility is completely smeared when both signal or idler beams are considered. On the 
other hand, we have demonstrated that for input Fock states, i.e. for sub-Poissonian 
light, our model can provide sub-Poissonian light governed by coupler parameters. 
However, for input coherent light, partially coherent light and chaotic light can be 
generated. This was demonstrated in terms of the second-order correlation function. 
Concerning the photon antibunching we have shown for both number and coherent 
input states that the outgoing field oscillates between bunching and antibunching giving 
a good indication that it need not be direct relation between sub-Poissonian statistics 
and photon antibunching as shown in the literature before |145II144] . 
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In this doctoral thesis we have studied the quantum properties of several models which 
have been classified as statical and dynamical systems. The first part has been de- 
voted to investigate the properties of the statical models including the superposition 
of squeezed displaced number states with and without thermal noise. Also we have 
developed a new type of multidimensional squeeze operator including two different 
squeezing mechanisms. In the second part the dynamical models were given to show 
the interaction between modes in the nonlinear optical coupler. The results of the 
statical models can be summarized as follows: 

1) We introduced general class of quantum states, as a single mode vibration of 
electromagnetic field, sudden squeezed-plus-displaced by a collection of two displace- 
ments 180*^ out of phase. For such class the quantum statistical quantities, such as 
the second-order correlation function, quasiprobability distribution functions and the 
photon-number distribution function are examined analytically and numerically. Moreover, 
the phase properties and amplitude squeezing and phase squeezing for such class have 
been examined in the framework of Pegg-Barnett formalism. Generalizations of some 
considerable results given in the literature earlier have been reported. Also we have 
suggested ways for generation of such superposition in the framework of micromaser 
and trapped ions. Regimes of strongly nonclassical behaviour of such states have been 
demonstrated. 

2) We studied the effect of thermal noise on the statistical behaviour of the super- 
position of squeezed and displaced number states. The main tool in this analysis is 
the s-parametrized characteristic function. For such a superpostion states an exact 
expressions for quasiprobability distribution functions, normalized second-order co- 
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relation function, photon-number distribution and phase space distribution have been 
obtained and discussed in detail. The behaviour of such system reavels that the ori- 
gin of the nonclassical effects is the correlation between different oscillators in phase 
space. 

3) We developed a new squeeze operator, which is related to the time-dependent evolu- 
tion operator for Hamiltonian representing mutual interaction between three different 
modes. Squeezing phenomenon as well as the variances of the photon-number sum 
and difference have been considered. Moreover, Glauber second-order correlation 
function, the quasiprobability distribution function and phase distribution for different 
states have been discussed. We have shown also that the origin of the nonclassical 
effects of this operator model is the correlation between the systems. 

On the other hand, the dynamical models were given in the second part to show 
the interaction between modes in the nonlinear optical coupler. We have given three 
types of such device and the results can be summarized as follows: 

4) We derived the quantum statistical and dynamical properties of nonlinear optical 
couplers composed of two nonlinear waveguides operating by the second subhar- 
monic generation, which are coupled linearly through evanescent waves and nonlin- 
early through nondegenerate optical parametric interaction. Main attention has been 
paid to generation and transmission of nonclassical light, based on a discussion of 
squeezing phenomenon and quasiprobability distribution functions. The beams have 
been initially considered in coherent and number states. In particular, results have 
been discussed in dependence on the strength of the nonlinear coupling relatively to 
the linear coupling. We have shown that if thermal fields enter initially the waveguides 
the coupler plays similar role as a microwave Josephson-junction parametric amplifier 
to generate squeezed thermal light. 

5) We showed that a nonlinear asymmetric directional coupler composed of a linear 
waveguide and a nonlinear waveguide operating by nondegenerate parametric ampli- 
fication is a source of single-mode squeezed light. This fact has been demonstrated, 
under certain conditions and for specific modes, for incident coherent beams in terms 
of the quasiprobability functions, photon-number distribution and phase distribution. 

6) We studied the quantum statistical properties of an optical field propagating in- 
side a directional coupler operating by nondegenerate parametric amplification. We 
have investigated the effect of switching between the input modes and the outgoing 
fields from the coupler. Particular attention has been paid to two mode squeezing 
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and second-order correlation function. Incident number and coherent states have 
been considered. Furthermore, regimes for generation and transmission of squeezed 
and/or sub-Poissonian light have been found. 
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Appendix A 



In this appendix we deduce the result of the Fourier transformation of equation ( |5.1.17| . 
This is basically depending on the generating function of the Laguerre polynomial. We 
restrict ourselves to the derivation of the result of one term for Wigner function, i.e. 
s = 0, however, the others can be found in a similar way. 

Now let us assume we have such type of integral 

In = I d^Cl^n{\kf)exp[-\k\^/2 + {C-C)a]exp(fiC -PX). (A.l) 

where the explicit form of k has been given in l |5.1.18| . Using the generating function 
of Laguerre polynomial |178| 



exp(-T^ 
1 - t 



(A.2) 



n=0 

in equation (A.l), we get 

oo 



^ nl 

1 



1 



1 - t 



1 - t 



d^Cexp 



d'^Cexp 
t + 1 



t + 1 



lfe|2 + C(a-j8*) + C*(j8-a) 



2(1-0' 

082cosh(2r) + (jS^ + j8*2)C^S^] + C(a - + C(fi -a)] ,(A.3) 



2(1 - t) 

where the notations and have the same meaning as before. Applying the fol 
lowing identity 1^ 

exp[-B|C|2 + (c/2)C*2 + {Ci/2)C^ + DiC + DC]d^C 
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K 



exp \ ^[DDiB + d2(ci/2) + D?(c/2) 



with 



K = - cci, ReK > 0, Re[B + (ci + c)/2] > 0, 
to (A.3), after minor algebra we arrive at 



E 

n=0 



T /" 



2jt 

TT7 



exp 



2:jt 
1 + t 



exp(-2z) exp 



, (1-0^ 
■ (1 + 

tz 



■(1 + f) 



2;Texp(-2z) ^-^Ln(4z 
. n! 



n=0 



(A.4) 
(A.5) 

(A.6a) 
(A.6b) 
(A.6c) 



(A.7) 



where 

z = |j8 - a|2cosh(2r) + - af + (fi* - af]SrCr, 

= (x - a)^exp(2r) + y^exp(-2r). 

In the transition from (A.6b) to (A.6c) we used the generating function of Laguerre 
polynomial (A.2) again, however, in (A.7) we used j8 = x + iy as before. Finally, (A.6c) 
shows that 

In = 2;T(-l)"exp(-2z)Ln(4z), (A.8) 

and this is the value of the required integral. For Q-function the right-hand side of 
(A.6) will be expressed in terms of the generating function of Hermite polynomials. 



Doctoral Thesis 



122 



Faisal Aly Aly El-Orany 



Appendix B 



In this appendix we derive the relations j5.3.8| . Before solving the problem it is con- 
venient to remind that the connection between Schrodinger's picture and Heisenberg's 
picture for any operator is 

Aj{t) = sHt)Aj{0)S{t), (B.l) 

where S(f) is the unitary operator of time development. So that the required relation 
can be found by solving Heisenberg's equations of motion. 

The Heisenberg equation of motion for any operator O is given by 

The equations of motion related to the effective Hamiltonian l |5.5.3t are 

^ = XiAl + X2AI (B.3a) 

^ = AiA[ + A3A3, (B.3b) 

^ = X2AI - A3A2. (B.3c) 

Equations (B.3) with their Hermitian conjugates give a close system of 6 differential 
equations with time-independent coefficients which can be solved, for instance, by 
means of Laplace transformation. Nevertheless, they can be solved straightforwardly 
by differentiting one of these equation two times and subsituting from the others. Let 
us restrict ourselves to (B.3a) which gives (after minor algebra) 

^-(A? + Ai-Ai)^^0, (B.4) 
and consequently the charactristic equation is 

y3 - (Af -h A| - A|)y = 0, (B.5) 



and hence the roots of (B.5) are y = 0,+z = +yAj + Ag - A§. Now considering the 
condition Af + Ag > A|, the solution of (B.4) is 

Ai{t) = ci + C2Cosh{tz) + C3sinh(fz), (B.6) 

where ci, C2, C3 are constants which can be specified from (B.6) and its derivative. After 
some simplification we arrive at 

Ai(f) = Ai(0)[cosh(fz) -h ^ sinh2(fz/2)] 



Doctoral Thesis 



123 



Faisal Aly Aly El-Orany 



Chapter 8. Appendices 



+AU0)[^smh{tz) - fl^sinh2(fz/2)] 
z 

+A3(0)[— sinh(fz) + sinh2(fz/2)l. (B.7) 

z 

Now using the substitution = fAy,; = 1,2,3, (B.7) reduces to the first equation in 
j5.3.8| where we wrote Aj{0) = Aj in these relations. Similarly, the relations related to 
Aj{t),i = 2,3 can be deduced. 
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In this appendix we give the explicit forms for the expectation values of cross 
photon-number operators between various modes for three-mode squeezed coherent 
states. The derivation is straightforward and the relations j5.3.8| should be frequently 
used: 

+ |(a2/2 + "3/3) ("2^1 + "3^2)1^ 
+ (jf\aif + ff - l)|a2gi + a3g2)|^ 

+ {\aif + 1) {\a2f + i)/2g3(gi/i + ^3/2) + +(|a3|^ + i)/3g3(g2/i + ^3/3) 

+hgz{f2gi\cx2f + /3g2|a3|^) + (aia2 + ata2)/iSf3(/igi + /2g3) 

+ (0203 + a2a^)gs{fihgi + fif2g2 + /2/3g3) + {ocicc^ + cclal)fig-5{fig2 + hg^) 

+fi\oc2g\ + a^g2f[{ocia2 + a^a2)/2 + (aias + ala^)h] 
+g-i\cc2f2 + a3/3|^[(aia2 + aXaDg^ + (aias + a\al)g2] 
+ (aia2 + a\al)h2[f if 2^2 + ffh-^ - h-^] 
+ (aia3 + a\al)h2\Jif2h2 + ffh^ - hi] 
+/ih2[(af + af)\a2\^f2h^ + (a? + af)\a-i\^hhi 

+ (afa2a3 + af a^a^)(hi/2 + h-^h)], (C.l) 

(n2n3) = h|g|[(|ai|2 + if - 3] + (|ai|2 + 1) g||hia3 + h3a2|^ 

+hl\gia2 + g2a3|^ + h2h-^gig-5{2\a2\^ + 1) 
+hih2g2g3(2|a3|2 + 1) + hzg^igihi + g-sh-s){aia2 + a^aa) 
+^2^3(^3^11 + g2^i2)(aia3 + oclal) 
+h2g-5{g2h-5 + gihi)(a2a3 + agas) 

+ |(gia2 + g2a3)(^iia3 + ^i3a2)|^ 

+g^\hia^ + h^a2f[gi{aia2 + ala^) + g2{ocia^ + octal)] 
+h2|gia2 + g2a5f[hi{aia^ + alal) + h3(aia2 + a^al)] 
+h2g^{hig-s + h2g2)(aia3 + oclal) 
+h2g^{h^g-5 + h2gi){aia2 + a^ag) 
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+h2g^{higi + h^g2){ala2a^ + afa^a^) 

+h2hzgigz{afal + afaf) + hih2gig^{afal + afaf). (C.2) 

Corresponding relation between modes 1 and 3 can be obtained from (C.l) using the 
transformation l |5.3.21| . 
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Appendix D 



In this appendix we give the derivation of the solution of equations of motion j6.1.2| . 
The Heisenberg equations of motion are 
ddi 



dt 
dt 



If t 

- = -i(i»idi - 2iAidj exp(if^ii) - iX-^dzexpl-icpiit)] - iA4d2 exp[i(/)2(f)], (D.al) 
= -i(L»2d2 - 2iA2d2exp(ifjLi2) - iAsdi exp[i(/)i(f)] - iA4d[exp[i(/)2(f)]. (D.a2) 



Substituting di = Aexp(^^ii) and d2 = Bexp(^^i2) in terms of the slowly varying forms 
A and B, equations (D.a) take the form 
dA 



dt 



dB 
It 



+ y)A - 2iAiA^ - iAsBexp 



-iA4B^exp 



-I K + ^4>2[tl 



-i{(^2 + - 2iA2B^ - iAsAexp 



.(;i2 -Mi)^ . , 
I n 



(D.bl) 



exp 

Now let us define 

Qi = A + Q2 = B + B^ Pi = A - p2 = B- BI 



(D.b2) 



(D.b3) 

Equations (D.bl)-(D.b2) cannot be solved directly and hence some restrictions should 
be considered, so that we shall consider 4>i{t) = 5(^2 - and 4>2{t) = ^{p.2 + P^i)t and 
thus these equations with their Hermitian conjugates lead to the following equations: 

dQi ..ft ..ft dQ2 



dt 
dp I 

~dr 

where 



-ik-Pi - iX-P2, 



-ik+Qi - iX+Q2, 



dt 

dp2 
dt 



= -i/_P2 - iA_Pi, 



-i/+Q2 - iA+Qi, 



1 1 

A± = A3 + A4, k+ = (A)i + -111 + 2Ai, /+ = (i»2 + -^112 ± 2A2. 



(D.cl) 



(D.c2) 



The second order differential equations for Qj can be obtained from (D.cl) after 
some minor algebra as 



where Qj 



+ QiQi = -yJgig2Q2, 
= y/gjQj, i = 1,2 and 



d^Q2 
df2 



+ Q.2Q2 = -VmrnQi, 



ip.di) 
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gi = k^X+ + A„/+, §2 = i?+A„ + A+/„, 

Qf = A_A+ + k^k+, Q'^ = A_A+ + /_/+. (D.d2) 

Now using the following transformations: 

Oi = XcosQ + VsinQ, O2 = VcosQ - XsinQ, (D./l) 

where X and V are two Hermitian operators, equations (D.dl) can be split into two 
equations as 



provided that 9 = ^tan"^ \ jj^f^^k+ )' where 



Qi = 
^2 = 



qI cos^ e + Ql sin^ - ^^1^2 sin(20) 



(D./2) 
(D./3) 

{D.fA) Equations (D.f2) can be 



cos^ + g|j^2 Q ygTgg sin(20) 

easily solved and consequently Qj, 7 = 1,2 can be obtained. Similar procedures can 
be done to get Pj, ; = 1,2 and eventually the expressions for dy(f), ; = 1,2 can be 
traced out in the form j6.1.2|l . The time dependent coefficients in l |6.1.2| are given by 



Ki(f) = Fi(f) - 



Mi(f) 

Ni{t) - 
whereas 



[k+ + fe-]Gi(f) + [A+^ + A_]S(f) 
9i 



92 



[fe+-fe_]Gi(f) + [A+^-A_]S(f) 
91 



1 

2 

1 

2 



1 + 



92 
9i 
92 
9i 



C{t) - i 
C{t) - i 



92 



K2(0 = F2(f) - - 



[A++A„]Gi(f) + [/+^+/_]S(f) 
9i 

9i 



92 



m) - 

M2(0 
N2(f) - 



[/++/_]G2(f) + [A+^+A_]S(f) 
9\ 



92, 



[/+ - /_]G2(f) + [A+ -A„^]S(f) 

9i 



1 

2 

1 

2 



1 + 



92 
9i 



C{t) - i 



92, 



92 
9i 



- 1 C(f) - i 



[A++A„]G2(f) + [fe++fe_^^]S(f) 

9i 

[A+ -A_]G2(f) + [fe+ -fe_^]S(f) 

^1 



(D.h) 



In the above equations we have defined 
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F2{t) = cos(ff^2) cos^ + cos(ff>i) sin^ 6. 



sm{tQi) 2 
Gi(n = z ^cos I 



+ — ^ — - sin^ I 

^2 



^ sinffQa) 9^ sin(fQi) . 9^ 
G2(f) = — cos^ + — sin^ 0, 



2 V ^2 



S(f) 



1 

2 



COs(ff22) 

sin(ff^2) 

f^2 



cos{tQi) sin(20), 
sin(ff^i) 



sm(20). 
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Appendix E 

In this appendix we give the exact expressions for the time dependent functions for 
the solutions of the Heisenberg equations of motion for interacting modes according 
to the Hamiltonian ( |6.3.1| . It is reasonable to mention that the solutions in this case 
can be deduced using the similar technique as that given in appendix D. Thus we have 

/i(0 = cos{tQi) cosh^ (f) - cosh(ff^2) sinh^ (p. (E.la) 

/2(0 = t[^-^^]sinh(2(/)) 

, , ,sin(ff^i) , 9 , sinh(ff^2) . , 9 ,n ,^ . ^ 

hW = Ai[ ^ " cosh^cj) -^-^sinh2(^], (E.lc) 

/4(f) = i[cos(ff2i) -cosh(ff^2)]sinh(2(/)), (E.ld) 

h,(f) = ^[!^-!i!^^]sinh(2e) 

, .sinhfff^/) , 9 „ sinfff^^) , 9 ,^ ^ v 

+Ai[ — ^cosh^e p-^sinh^Q], (E.2a) 

1 

h2(f) = -[cosh(ff24) - cos(fQ3)]sinh(20), (E.2b) 

h^{t) = cosh(fQ4) cosh^ - cos{m^) sinh^ 6, (E.2c) 

A3 sin(ff23)sinh(ff24) . 
^4(0 = y[ — ]sinh(20), (E.2d) 

where 

Qi = [kf cosh^ (/> + A| sinh^ - A2A3 sinh(2(/))] ^ , (E.3a) 

^^2 = [kf sinh^ + A§ cosh^ (p - A2A3 sinh(2(/))] ^ , (E.3b) 
and 

^-^tanh-^ U'y\2 )> 



with fei = y Ag - A^. 

The other coefficients can be obtained with the help of the substitution Ai < — > A3 
and hence 

Qi < — > Q-^, Q2 ^ — > 4* ^ — ^ 

{h{i).f2{i).h{i).h{i)) ^ {g2{i).gi{i).gdi).gm. 

{h,{t),h2{t).h-i{t).hm ^ {W)M{t).m).h{t))- (E.5) 
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Sf rucne shrnuf i vysledku diserf acni prace (summarsz in Czech) 

V teto disertacni praci jsou studovany kvantove vlastnosti nekterych statickych a dy- 
namickych systemu. Prvm' cast je venovana vyzkumu vlastnosti' statickych modelu, 
zahrnuji'cim superpozice stlacenych a posunutych Fockovych stavu s termalm'ch sumem 
a bez termalniho sumu. Je rovnez zaveden novy typ vi'cerozmerneho operatoru stlaceni, 
ktery zahrnuje dva ruzne typy mechanizmu stlaceni. V druhe casti jsou vysetrovany 
dynamicke modely v souvislosti s interakci modu v nelinearni'ch optickych vazebnich 
prvcich. Motivace je dana dfi'vejsim vyzkumem techto prvku jako slibnych zanzem 
pro generaci neklasickeho svetla. Jsou diskutovany tri typy techto zanzem, zalozene 
na interakci dvou, tri a ctyf modu pf i sifeni prostf edim. Jsou urceny rezimy generace 
a pfenosu neklasickeho svetla. 
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